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A super-operator  formulation  for  decoupling  the  symmetry  blocks  of 
the  full  2p-h  TDA  electron  propagator  matrix  is  presented.  The  accuracy 
of  this  decoupling  is  explored  by  using  the  full  2p-h  TDA  overlap 
amplitudes  to  calculate  relative  ionization  intensities  for  N2  and  H^O. 
The  effect  of  excluding  orbitals  on  the  calculated  ionization  potentials 
and  relative  ionization  cross-sections  is  reported.  The  use  of  intensi- 
ties obtained  from  successively  larger  subsets  for  labeling  the  photo- 
electron spectrum  is  suggested. 

Expressions  are  formulated  for  the  correlated  treatment  of  photoion- 
ization cross-sections  using  Feynman-Dyson  amplitudes  of  the  electron 
propagator  to  describe  the  ionization  processes  and  an  orthogonal i zed 
plane  wave  description  for  the  photoelectron.  Explicit  equations  are 
derived  for  the  theoretical  treatment  of  experimental  cross-sections 
using  either  spherical  grid  deflection  type  spectrometers  or  those  using 
a particular  collection  angle.  A scheme  is  suggested  for  systematic 
improvement  of  the  theoretical  description. 
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A dilated  electron  propagator  for  the  direct  calculation  of  reso- 
nance energies  and  widths  for  quasi -bound  anions  of  closed  shell  atoms 
is  proposed.  The  treatment  employs  creation  and  annihilation  operators 
defined  on  bi-variationally  obtained  bi-orthogonal  orbital  basis  set. 
Explicit  equations  for  the  dilated  electron  propagator  correct  through 
second  order  in  electron  interaction  are  presented.  Higher  order  and 
renormalized  treatments  are  suggested  and  the  details  of  a computational 
strategy  outlined.  The  bi -variational  SCF  scheme  for  di lata ted  closed 
shell  atomic  Hamiltonians  is  studied  and  discussed.  Applications  to 
the  Be  atom  are  presented. 
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INTRODUCTION 


The  bound-continuum  transitions  which  involve  the  ejection  of  elec- 
trons from  bound  states  into  continuum  or  the  attachment  of  electrons  in 
the  continuum  to  bound  or  quasi-bound  states,  due  to  the  quantal  nature 
of  their  energetics,  provide  a definitive  insight  into  the  electronic 
structure  and  dynamics  of  atoms  and  molecules.  Since  the  attachment  or 
detachment  of  electrons  may  lead  to  substantial  realignment  of  inter- 
electronic  interactions,  the  independent  particle  schemes  which  take 
account  of  these  interactions  only  in  a particular  average  way  are  not 
always  satisfactory.  As  such,  correlated  treatments  of  bound-continuum 
transitions  which  allow  for  dynamic  correlation  of  interelectronic  inter- 
actions are  of  paramount  importance  in  attempting  a theoretical  descrip- 
tion of  chemical  phenomena. 

In  its  spectral  representation,  poles  of  the  electron  propagator 
(Linderberg  and  Ohrn,  1973) 


(E) 


lim  Y 
n-HD  s 


fs(p)fj(q) 

E-{ES ( N+l ) -Eq (N) }+in 


+ 


gs(p)g|(q) 

E— (Eq(N) — E$ (N— I ) }— in 


where  Feynman-Dyson  overlap  amplitudes  f$  and  g$  are  defined  as 

fs(p)  = <N|a  |N+l,s> 

gjp)  = <N-1  ,s | a | N> 
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permit  the  direct  calculation  of  transition  energies  involved  in 
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attachment  {Es (N+l ) -EQ ( N ) } and  detachment  {EQ(N)-E  (N-l)}  of  electrons  in 
atoms  and  molecules.  The  dynamic  correlation  of  inter-electronic  inter- 
actions is  incorporated  through  summation  over  all  possible  (N+l,s)  or 
(N-l,s)  electron  states  that  can  emerge  from  the  initial  ground  state  of 
energy  EQ(N).  Other  correlated  treatments  of  bound-continuum  transitions 
necessitate  the  calculation  of  total  energy  for  both  ground  (Eq (N ) } and 
the  final  state  { E ( N+l ) or  E(N-l)},  and  transition  energy  is  calculated 
as  the  difference.  These  latter  methods  suffer  from  the  significant  loss 
of  accuracy  inherent  in  the  subtraction  of  two  nearly  equal  total  energies 
to  obtain  a much  smaller  transition  energy.  Also,  the  approximations 
resorted  to  for  computational  tractability  of  these  methods  do  not  always 
lend  themselves  to  cogent  mathematical  formulations.  Care  must  also  be 
taken  to  avoid  disparate  levels  of  approximating  electron  correlation  in 
the  calculation  of  each  state.  The  electron  propagator  formalism  does 
not  suffer  from  these  defects,  and  various  approximations  lend  themselves 
to  systematic  improvement.  This  has  prompted  us  to  formulate  theoretical 
schemes  which  should  make  possible  an  accurate  treatment  of  bound-contin- 
uum transitions  using  the  electron  propagator  formalism  in  conjunction 
with  the  ideas  of  dilated  hamiltonians  (Aguilar  and  Combes,  1971;  Balslev 
and  Combes,  1971;  Simon,  1972). 

The  causal  electron  propagator  (Linderberg  and  Ohrn,  1973;  Simons, 

1976) 

G(x,x',x)  = -ie(x)<iKx,x)^+(x' ,0)>  + i0(-x)<^+(x'  ,0)iMx,t)> 

where  0(x)  is  the  Heaviside  step  function,  and  are  the  electron 
field  operators  in  the  Heisenberg  representation.  The  brackets  define  a 
ground  state  average  relative  to  some  normalized  density  operator  p 
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<^(x,x)4^(x 1 ,0)>  = Tr{pi|>(x,x)t|/*"(x 1 ,0)} 


Introducing  a discrete,  orthonormal  set  of  spin  orbitals  (u^x)}, 
the  electron  field  operators  can  be  expanded  in  a set  of  spin  orbital 
annihilation  and  creation  operators. 


asM 


u*(x)^(x,t)dx 


In  this  representation  the  electron  propagator  is  expressed  as  a 
matrix  with  the  elements 

Gsr(x)  = -ie(T)<as(x)aJ(0)>  + i0(-x)<aJ(O)as(x)> 


The  Fourier  transform  of  the  electron  propagator  satisfies  the  following 
equation  of  motion: 


(E) 


E'1<[as,a];]+>  + E‘2<[[as,H]_,a+]+> 
E'3<[[[as,H]_,H]_,aJ]+>  + ... 


where  H is  the  Hamiltonian 


H = 


hi 


<rr ' 


rr 1 ,ss 1 


4.  4. 

I ss 1 >a  a ,a  , a 
rrs  s 


and 

<rr ' | | ss  * > = [u*(l)u*, (2)r-1 (1-P  )u  (l)u (2)dxdx 

) ' ‘ 12  12  5 S 12 

By  introducing  a linear  space  L = (a^ , a ^a^am  ( item ) , ...}  of  fermion 
like  field  operators  which  supports  a scalar  product 

(Xi  IXj)  = Tr{p[xt,Xj]+} , Xi,XjeL 
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by  defining  the  identity  superoperator  i and  the  Hamiltonian  superoperator 
H as 


IX  = X 

HX  = [X,H]_  X e L 

and,  finally,  by  arranging  the  annihilation  operators  in  a row  vector  and 
the  creation  operators  in  a column  vector,  the  equation  of  motion  for 
the  Fourier  transformed  causal  electron  propagator  may  be  recast  as 
(Goscinski  and  Lukman,  1970;  Pickup  and  Goscinski,  1973) 

G(E)  = ( a | (EI-H) _1 a) 

which  may  be  approximated  via  the  inner-projection  technique  (Lowdin, 
1965)  as 

G(E)  - (a|h)(h|(EI-H)h)_1(h|a) 

where  the  choice  of  projection  manifold  h defines  a particular  approxima- 
tion to  the  electron  propagator. 

In  chapter  one  we  describe  the  use  of  symmetry  adapted  linear 
combination  of  annihilation  and  creation  operator  products  to  decouple 
the  symmetry  blocks  of  the  electron  propagator  matrix.  We  demonstrate 
its  use  for  making  manageable  the  computation  of  transition  energies  for 
ionization  processes  using  the  Tamm-Dancoff  approximation  (TDA) . These 
results  are  used  for  comparative  investigation  of  TDA  and  other  approxi- 
mants  (Mishra  and  Ohrn,  1980a). 

The  necessary  formalism  for  the  use  of  Feynman-Dyson  amplitudes  of 
the  electron  propagator  for  correlated  treatment  of  photoionization 
cross-sections  is  fabricated  in  chapter  two  (Mishra  and  Ohrn,  1980b). 
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The  results  of  cross-section  calculations  using  this  formalism  with  plane 
wave  approximation  for  the  photoelectron  are  used  for  further  examination 
of  the  TDA  propagator. 

In  chapter  three,  the  electron  propagator  equations  are  generalized 
to  permit  the  description  of  transitions  involving  quasi -bound  states 
described  by  non-hermitian  Hamiltonians.  The  use  of  this  more  general 
dilated  electron  propagator  (Mishra,  Froelich  and  Ohrn,  1980)  for  a 
unified  treatment  of  all  bound-continuum  transitions  involving  bound 
states  or  quasi-bound  states  for  the  calculation  of  ionization  energy, 
electron  affinity,  shape,  and  Feshbach  resonances  (Taylor,  1970)  for  atoms 
and  molecules  is  adumbrated,  and  the  results  for  zeroth  order  dilated 
electron  propagator  (Mishra,  Ohrn  and  Froelich,  1980)  are  analyzed  and 
discussed  in  this  chapter. 


CHAPTER  1 


BOUND  CONTINUUM  TRANSITIONS  WITHIN  2p-h  TAMM-DANCOFF  APPROXIMATION 


1.1  The  2p-h  Tamm-Dancoff  Approximation,  (2p-h  TDA) 

We  have  noted  earlier  that,  in  the  superoperator  formulation 
(Goscinski  and  Lukman,  1970;  Pickup  and  Goscinski,  1973),  the  electron 
propagator  matrix  is  given  by 

G(E)  = (a | (EI-H) “ 'a)  (1.1) 

which  may  be  approximated  via  the  inner  projection  technique  (Lowdin, 
1965) 


G(E)  ~ (a|h)(h| ( El -H ) h ) " 1 ( h | a ) 


(1.2) 


It  has  been  shown  (Born  and  Ohrn,  1979)  that  the  choice  of  projection 
manifold  h as 

h = hj  © h3  (1.3) 

where 


{a.} 


-3 


{W 


k£m 

^k£m 


= Nkfa,iakaA  + 5km<nk>ail  ' 'Wam> 
= nk{(l-<nt>-<nm>)  + <nJ,xnm>} 


<n^>  = occupation  number  for  the  spin  orbital  k 
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leads  to  what  has  been  labelled  as  the  2p-h  (2  particle-hole)  Tamm- 
Dancoff  Approximation  (Schirmer  and  Cederbaum,  1978a). 

With  this  choice  of  the  projection  manifold,  the  electron  propagator 
may  be  written  as 

G"‘(E)  - (hjlEI-Hjhj)  - (h1|Hh3)(h3|(E!-H)h3)'I(h3|Hh1)  (1.4) 

= - i (1.5) 

where 

2 = (h1|Hh3)(h3|(EI-H)h3)'1(h3|Hh1)  (1.6) 

is  the  so-called  self-energy  matrix  (Dyson,  1949). 

The  computational  bottleneck  in  implementing  the  full  2p-h  TDA  and 
other  self-energy  approximants  has  been  the  size  of  the  hk£m  manifold 
(-M3,  where  M is  the  size  of  the  basis)  and  the  consequent  enormity  of 
the  matrix  inversion  problem  in  eq.  (1.6).  Previously,  this  has  neces- 
sitated the  approximation  of  the  inverse  matrix  in  (1.6)  by  its  diagonal 
elements  (the  diagonal  2p-h  TDA)  or  the  approximation  of  H in  the  inverse 

A 

matrix  by  H0  (the  second  order  approximation).  These  approximations 
obliterate  correlation  effects  which  provide  a deeper  insight  into  the 
electronic  structure  of  chemical  systems,  e.g.  a breakdown  of  the  orbital 
picture  of  ionization  for  the  inner  valence  orbitals  (Schirmer  and 
Cederbaum,  1978b;  Mishra  and  Ohrn,  1980a). 

In  the  next  section,  we  describe  the  use  of  symmetry  adapted  linear 
combination  of  the  operator  products  in  the  h3  manifold  to  reduce  the 
size  of  the  matrix  inversion  problem  by  decoupling  the  evaluation  of 
electron  propagator  matrix  into  those  of  its  smaller  symmetry  blocks. 

This  makes  possible  the  evaluation  of  the  full  2p-h  TDA  propagator. 
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While  the  explicit  use  of  this  technique  has  been  made  only  for  the 
reduction  of  the  h^m  manifold  into  its  invariant  subspaces,  the 


significant  savings  in  the  higher  order  decouplings  of  the  electron 
propagator  (Redmon,  Purvis  and  Ohrn,  1975). 

1.2  Decoupling  of  the  Symmetry  Blocks  of  the  2p-h  TDA  Propagator 

The  reduction  of  h^^  manifold  may  be  achieved  by  using  the  decompo- 
sition 


where  the  summation  is  over  all  the  irreducible  representations  of  the 


ksym,  Jlsym,  and  msym  being  symmetry  labels  for  the  orbitals  k,  % and  m 
obtained  from  the  SCF  procedure.  We  limit  ourselves  to  those  point  groups 
which  have  only  one  dimensional  irreducible  representations . 


I © h\ 
i J 


(1.7) 


i = ksym  ® £sym  ® msym 


(1.8) 


Obviously 


(1.9) 


This  affords  the  rewriting  of  E(E)  as 


2(E)  = l © Zn(E) 


(1.10) 


where 


51 (E)  = (h{|Hhj)(hj|(EI-H)h^)‘*(hj|flhj) 


(1.11) 


9 


We  further  reduce  the  size  of  the  problem  by  realizing  that  the 
creation  and  annihilation  operators  are  tensor  operators  of  rank  h with 
respect  to  the  spin  angular  momentum,  and  as  such,  the  following  decompo- 
sition 


hg  = q1  ® f1  © d1 


(1.12) 


may  be  achieved  by  using  the  Wigner-Racah  Tensor  calculus  (Linderberg  and 
Ohrn,  1973;  Messiah,  1958)  whereby 


^akK^aUamy^Ms^Ms 


(1.13) 


= l (-)Ms'(2s'+irS(2s+l  ^(Xp-M1.) 

<»Ms,,X,u  s 

v ,h  s'  s , + - - 

K Ms'  "Ms)akKa£Xamy 


V = <->*Vy 
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^qk£m^ 3 
^qkilm^V2 


% = 


^aki<aJUamy)  3/2 


^aki<aj!,Aamy  ^ 1/z 


M~%  + 

K&na  kaail8am8 


k£m r + + 

/2  1 kB  &B  mB  kcta£Bama 


(1.14) 

(1.15) 


akaa£aamB^ 


(1.16) 


( q ) V2 


(akAxSmp>-t 


N"'5 

k£m  r + + 

/3  ^akaa£aama  akBa£aamB 


“ak6a£BamcJ 


(1.17) 


3/2 


^ ^ L'  0 m ^ — 3/ 
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-}2  + 
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(1.18) 
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k&m  r + + . 

/2  lakaail6ama‘akaatoam6"6k£<nk>am6 

-5km<nk>alB}  l<m  (1'19) 

NktonfakaalBato''5kJ<nk>al!.e^  d-213) 

l = m 

_u 

n z 

ktof,  + + 

~jz  l<:akBa«BamB  akaa«.Bama 

+akaataamB-36kH<nk>amB 
+3Skni<nk>atB1  l<m  (1'21) 

In  further  discussions  we  drop  the  superscript  for  the  irrep  of  the 
point  group,  since  the  reduction  of  the  projection  manifold  has  uncoupled 
the  different  symmetry  blocks,  and  we  treat  them  separately.  Thus,  for 
every  symmetry  block 


E(E)  = [(a|vq)(a|vf)(a|vd)] 

(d I (EI-H)q) (q | (El-fi)f ) (q [ (EI-H)d) 

-1 

(q 1 va) 

(f | (EI-H)q) (f | (EI-H)f) (f | (EI-H)d) 

(f | va) 

(d | (EI-H)q) (d | (EI-H)f ) (d | (Et-H)d) 

(d | va) 

(1.22) 


and  since 

(a|vq),  (q | va) , (f | (EI-H)q) , (d J (Ef-fi)q) , (q| (EI-H)f ) , 
(q | (Ef-fi)d) 

are  null  matrices, 


E(E)  = C(a|vf)(a|vd)] 

(f | (EI-H)f)(f | (EI-H) |d) 

-1 

(f |va) 

(d  (EI-H)f ) (d  (EI-H)  d) 

(d  va) 

k£m 


(ak<(atXamu>S)i 


k£m 


- (4<a)U5inp>Msi^ 


(1.24) 
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Diagonal ization  of  the  blocked  inverse  matrix  of  Eq.  1.24  leads  to 


Eij(E) 


- i 

n 


m.(n)mj(n) 
~ 1 E=iT 


(1.25) 


Detailed  expressions  for  mn.  (n)  and  in  terms  of  two-electron 
integrals  and  orbital  energies  may  easily  be  obtained  using  relations 
provided  by  Born  and  Ohrn  (1979)  employing  the  HF  density  operator  and 
integral  occupation  numbers.  These  are  identical  to  those  of  Schirmer 
and  Cederbaum  (1978a). 


1.3  Comparison  of  the  Full  2p-h  TDA  and  Other  Self-energy  Approximants 
As  introduced  earlier,  the  electron  propagator  in  its  spectral 
representation  is  given  as  (Linderberg  and  Ohrn,  1973). 


(E) 


1 im  l 

rpKH  s 


fs(p)f*(q) 


E+E 


0(N)-Es(N+l)+in 


gs(p)g*(q)  i 
E+Es(N-l)-E0(N)-inJ 


(1.26) 


The  quantities  needed  for  the  comparison  of  different  decouplings  are  the 
ionization  energies  and  the  Feynman-Dyson  overlap  amplitudes  fs  and  gs 
defined  as 


fs(p)  = <N|ap|N+l,s>  (1.27) 

9S(P)  = <N- 1 , s | ap | N>  (1.28) 

Schirmer  and  Cederbaum  (1978b)  have  used  the  full  2p-h  Taimi-Dancoff 
Approximation  to  the  electron  propagator  to  investigate  the  inner  valence 
ionization  of  several  molecules  and  have  compared  its  accuracy  with  the 
diagonal  approximations.  An  interesting  feature  that  has  emerged  from 
their  investigations  is  the  complete  breakdown  of  the  orbital  picture  of 
ionization  for  the  inner  valence  orbitals.  However,  their  comparative 


12 


investigation  of  the  full  2p-h  TDA  has  been  done  using  only  a subset  of 
the  one  particle  orbitals  (Schirmer  and  Cederbaum,  1978b;  Cederbaum  et 
al . , 1978)  and  by  comparing  the  poles  of  the  propagator  with  experimental 
and  other  theoretical  values  for  the  ionization  energies.  This  criterion 
represents  a particular  bias,  since  it  does  not  reflect  the  quality  of 
the  overlap  amplitudes.  The  Feynman-Dyson  amplitudes  determine  the 
spectral  density  function  (Linderberg  and  Ohrn,  1973), 


A(x,x‘ ;E) 


Tfk(x)fJ(x')6(E-Ek) 

k 

Igk(x)g*(x')6(E-Ek) 

'■  k 


E > y 

E < y 


(1.29) 


where  the  parameter  y in  this  case  is  the  chemical  potential  of  the 
N-electron  ground  state.  Since  a good  deal  of  important  information  can 
be  extracted  from  the  spectral  density  function,  e.g.  the  one-particle 
reduced  density  matrix  (Linderberg  and  Ohrn,  1973), 


y(x,x  1 ) 


c A(x,x ' ;E)dE 


(1.30) 


it  is  important  to  choose  a decoupling  approximation  which  not  only  yields 
accurate  ionization  energies  but  also  an  accurate  spectral  density.  In 
this  section  we  delineate  a scheme  for  obtaining  the  Feynman-Dyson 
amplitudes  and  ionization  energies  from  the  2p-h  TDA  propagator.  The 
ionization  energies  so  obtained  are  compared  to  those  obtained  from  other 
approximations.  In  the  next  chapter,  we  present  expressions  relating  the 
overlap  amplitudes  to  ionization  cross-sections.  The  calculated  photo- 
ionization intensities  will  then  be  used  for  further  comparison  of  the 
full  2p-h  TDA  with  other  approximations. 
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To  extract  the  Feynman-Dyson  amplitudes,  following  Purvis  and  Ohrn 
(1974)  we  recast  the  Dyson  equation  (1.5)  as  the  operator  equation 


G-1  (E)  = E - L(E) 


(1.31) 


or  equivalently, 

(E-L(E) )G(E)  = I (1.32) 

A formal  solution  of  this  equation  can  be  expressed  in  terms  of  the 
biorthogonal  expansion  (Layzer,  1963)  as 


G(E) 


$n(E)$*(E) 

E-'VE> 


(1.33) 


where  $n  and  Wn(E)  are  the  solutions  of  the  eigenvalue  equation 


L(E)$n(E)  = Wn(E)$n(E) 


(1.34) 


The  solution  of  this  equation  in  a discrete  basis  may  be  related  to  the 
overlap  amplitudes  from  elementary  residue  calculus  as 


lim  (E-EJG(E)  = 
E->E_  s 


Vs 


= r t (E  W*(E  ) 
S S^  s'  s'  s' 


(1.35) 


with 


,-i 


dW$(E)" 

dE- 


E=E. 


(1.36) 


and  £s's  denote  either  the  f or  g amplitudes  (Eq.  1.27  and  1.28).  A 
comparison  of  Eq.  1.26  and  1.33  reveals  that  the  ionization  energies  are 
solutions  of  the  equation 


Wn(E) 


E 


(1.37) 
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Table  I.  Full  2p-h  TDA  Ionization  Energies  for  N2 


Orbital 

ib 

Ic 

.full 

sdiag. 

I .P. (exp. ^ 

43.07 

43.0 

44.50 

41.90 

41.78 

41.77 

2CTg 

41.08 

40.76 

40.76 

38.89 

38.58 

38.39 

38.38 

39.3 

37.3 

29.25 

29.23 

29.17 

29.17 

37.13 

36.84 

36.84 

2au 

24.70 

24.61 

24.59 

24.59 

17.22 

17.11 

17.13 

17.13 

18.6 

18.6 

l7ru 

17.06 

17.26 

17.21 

17.21 

17.4 

16.8 

3ag 

13.92 

14.25 

14.2 

14.2 

15.1 

15.5 

a The  two  deepest  core  and  the  four  highest  virtuals  with  orbital 
energies  differing  by  1.0  H and  more  from  the  valence  region 
(-1.70  H - 0.0  H)  were  dropped  in  the  solution  of  Eq.  1.24. 


b The  two  deepest  core  and  the  three  highest  virtuals  with  orbital 
energies  differing  by  1.50  H and  more  from  the  valence  region  were 
dropped. 

c The  two  deepest  core  and  the  two  highest  virtuals  with  orbital 
energies  differing  by  2.0  H and  more  from  the  valence  reqion  were 
dropped. 


d Purvis  and  Ohrn,  1975. 
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The  equations  1.24,  1.35  and  1.37  were  solved  for  N2,  H20  and  CN"  by 
modifying  the  GRNFNC  program  (Purvis,  1973;  Born,  1979).  One-  and  two- 
electron  integrals  over  the  contracted  gaussian  type  orbitals  were 
generated  using  MOLECULE  (Almlof,  1974). 

The  calculated  and  experimental  ionization  energies  using  Dunning's 
(1970)  [9s,  5p/4s,  2p]  basis  set  for  N2  are  presented  in  Table  I.  This 
basis  set  is  comparable  to  that  employed  by  Schirmer  et  al . (1977),  and 
our  results  do  not  differ  in  any  significant  manner.  Surprisingly,  while 
the  diagonal  approximation  gives  the  correct  ordering  of  2 au  and  1ttu 
ionization  energies,  the  full  2p-h  TDA  ordering  is  marginally  in  error. 

The  full  2p-h  TDA  calculations  reveal  a complete  breakdown  of  the  orbital 
picture  for  the  2a^  orbital.  In  calculations  with  various  subsets  of 
orbitals,  we  find  that  the  smallest  subset  gives  results  which  correct 
the  wrong  ordering  of  2 a^,  1ttu  ionization  energies.  This  uncertain 
quality  of  getting  plausible  results  in  good  agreement  with  the  experimen- 
tal values  with  a small  subset  of  molecular  orbitals  underscores  the  need 
for  performing  such  calculations  with  successively  larger  subsets  until  a 
saturation  is  apparent.  In  the  present  calculation,  such  a saturation 
was  achieved  by  including  all  the  orbitals  with  orbital  energies  within 
2.0  H of  the  region  of  interest  (-1.70  - 0.0  H). 

The  results  from  the  calculations  on  the  H20  molecule  with  14  CGTO's 
and  26  CGTO's  are  presented  in  Tables  II  and  III.  The  14  CGT0  basis  set 
(Kurtz  and  Ohrn,  1979)  is  found  to  lack  flexibility  and  obliterates  the 
correlation  effects  which  should  lead  to  the  breakdown  of  the  orbital 
picture  of  ionization  for  the  inner  valence  2a  ^ orbital  as  shown  by  the 
calculations  with  the  26  CGT0  basis  set  (Kurtz  and  Ohrn,  1978).  The 
basis  set  dependence  is  typified  by  the  lbj.  principal  pole's  shifting 
nearly  2 electron  volts,  from  13.1  ev  with  14  CGTO's  to  11.3  ev  with 
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Table  II. 

Full  2p-h  TDA 

Ionization  Energies 

for  H20  (14  CGTO) 

Orbital 

2full 

zdiag. 

I .P. (exp. )a 

2ai 

36.4 

34.3 

32.2 

32.7 

3a^ 

14.8 

13.2 

14.7 

lb2 

18.9 

18.2 

18.6 

lbi 

13.1 

11.2 

12.6 

a Siegbahn  et  al . (1969). 
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Table  III.  Full  2p-h  TDA  Ionization  Energies  for  F^O  (26  CGTO) 


Orbital 

za 

Eb 

ZC 

£d 

zfull 

zdiag. 

I .P. (exp. )e 

38.12 

38.02 

37.88 

37.7 

2ai 

34.07 

34.05 

33.96 

33.8 

33.98 

33.08 

32.83 

32.54 

32.4 

33.4 

32.2 

32.57 

30.64 

30.60 

30.28 

30.20 

3a^ 

13.69 

13.45 

13.45 

13.49 

13.6 

13.7 

14.7 

lb2 

17.67 

17.89 

17.92 

18.03 

18.1 

18.2 

18.6 

lbi 

11.51 

11.12 

11.15 

11.25 

11.3 

11.4 

12.6 

a The  deepest  core  and  sixteen  highest  virturals  with  orbital  energies 
differing  by  1.0  H and  more  from  the  valence  region  (-1.50  - 0.0  H) 
were  dropped. 


b The  deepest  core  and  thirteen  highest  virtual s with  orbital  energies 
differing  by  1.50  H and  more  from  the  valence  region  were  dropped. 

c The  deepest  core  and  eleven  highest  virtual s with  orbital  energies 
differing  by  2.0  H and  more  from  the  valence  region  were  dropped. 

d The  deepest  core  and  six  highest  virtual s with  orbital  energies 
differing  by  3.0  H and  more  from  the  valence  region  were  dropped. 

e Siegbahn  et  al . (1969) . 
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26  CGTO's,  while  the  corresponding  change  in  the  diagonal  2p-h  TDA  is 
only  0.2  eV.  Once  again,  the  results  obtained  by  including  only  those 
orbitals  with  orbital  energies  within  1.0  H of  the  valence  region 
(-1.05  H - 0.0  H)  are  misleading.  Many  calculations  (Schirmer  and 
Cederbaum,  1978b)  have  appeared  in  the  literature,  where  all  orbitals 
differing  by  more  than  0.50  H from  the  valence  region  were  dropped  for 
convenience  from  the  calculations.  Our  results  indicate  the  need  for 
further  justification  of  such  procedures. 

Recently,  Ortiz  and  Ohrn  (1979)  have  used  the  self-energy  approximant 
of  Kurtz  and  Ohrn  (1979)  to  calculate  the  electron  affinity  of  CN"  as  the 
ionization  potential  of  the  5a  orbital.  Their  results  indicate  that  the 
description  of  many-body  effects  is  markedly  sensitive  to  the  quality  of 
the  computational  basis  set  employed.  To  understand  the  effect  of 
excluding  orbitals  obtained  from  large  computational  basis  sets  on  the 
full  2p-h  TDA  calculations,  we  present  in  Table  IV,  the  results  obtained 
for  the  5a  orbital  by  using  their  basis  sets  and  excluding  orbitals  in 
the  solution  of  Eq.  1.24.  The  fact  that  the  smaller  subset  of  the 
46  CGT0  M0 1 s gives  5a  ionization  potential  in  better  agreement  with  the 
experimental  value  illustrates  the  care  that  one  must  take  in  quantitative 
comparisons  when  limited  orbital  basis  sets  are  used. 

Since  ionizations  from  inner  valence  orbitals  involve  pronounced 
relaxation  and  intricate  correlation,  it  should  be  no  surprise  that  the 
diagonal  approximations  are  unable  to  unveil  a complex  process  like  break- 
down of  the  orbital  picture  of  ionization  for  inner  valence  orbitals. 
Non-diagonal  approximations  like  the  full  2p-h  TDA  provide  for  the  mixing 
of  many  more  configurations  and  should  be  more  efficacious.  The  results 
obtained  confirm  this  hypothesis  and  caution  against  the  spurious  results 
which  may  surface  if  the  choice  of  orbital  basis  set  is  not  balanced 
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Further  saturation  of  the  orbital  basis  set  was  not  attempted,  since  inclusion  of  more  than  twenty-fi 
orbitals  in  the  full  2p-h  TDA  calculations  is  not  feasible  with  the  present  facilities. 
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appropriately.  The  super-operator  formulation  for  the  use  of  symmetry 
adapted  linear  combination  of  creation  and  annihilation  operators  to 
block-diagonal ize  the  super-operator  resolvent  matrix  promises  considerable 
savings  in  the  implementation  of  all  propagator  approximations.  The  use 
of  successively  larger  subsets  in  the  full  2p-h  TDA  calculations  not 
only  provides  a criterion  for  the  reliability  of  these  results  but  reveals 
a clear  pattern  in  the  apportioning  of  intensities  between  the  different 
lines  into  which  the  principal  peak  of  the  orbital  picture  is  split  and 
as  such  provides  a sound  scheme  for  labelling  the  peaks  of  the  photoelec- 
tron spectrum.  The  necessary  formalism  and  calculations  are  described  in 
the  next  chapter. 


CHAPTER  2 


USE  OF  PROPAGATOR  AMPLITUDES  IN  THE  CALCULATION 
OF  PHOTOIONIZATION  CROSS-SECTIONS 


2.1  Introduction 

The  role  played  by  Einstein's  interpretation  of  photoelectric  effect 
in  the  understanding  of  atomic  and  molecular  structure  is  well  known. 
Since  then,  photoelectron  spectroscopy,  which  involves  an  analysis  of  the 
kinetic  energy  of  electrons  ejected  from  a system  by  an  incident  mono- 
chromatic beam  of  photons,  has  provided  most  suggestive  and  convincing 
evidence  for  the  orbital  picture  of  atoms  and  molecules  (and  also  its 
breakdown  in  the  most  naive  form).  A typical  photoelectron  spectrum 
consists  of  a series  of  bands.  If  an  orbital  description  of  the  elec- 
tronic structure  of  atoms  and  molecules  is  assumed,  the  first  band, 
corresponding  to  the  highest  observed  kinetic  energy,  is  associated  with 
ionization  from  highest  occupied  orbital,  having  the  lowest  ionization 
potential.  Successive  bands  then  correspond  to  ionizations  from  succes- 
sively lower  lying  orbitals,  the  spectrum  terminating  when  there  are  no 
further  levels  from  which  electrons  can  be  ejected  by  the  fixed  energy 
incident  photons.  With  the  advent  of  high  resolution  electrostatic 
kinetic  energy  analyzers  the  photoelectron  spectrum  reveals  a rich 
structure  of  satellite  peaks.  These  additional  peaks  arise  primarily 
from  "two-electron  processes"  involving  the  ionization  of  one  electron 
accompanied  by  excitation  of  a second  electron  and  provide  important 
information  about  electronic  correlation. 
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A scheme  for  understanding  the  involved  energetics  was  presented  in 
the  previous  chapter.  However,  the  real  power  of  photoelectron  spectros- 
copy as  an  analytical  tool  for  the  investigation  of  atomic  and  molecular 
structure  lies  in  the  details  of  the  principal  and  shake-up  spectral 
bands.  The  shapes  of  these  bands  are  related  to  photoionization  cross- 
sections  and  serve  as  fingerprints  for  atoms  and  molecules.  Comparison 
of  experimental  photoionization  intensities  with  calculated  cross-sections 
constitutes  an  important  assignment  aid  in  photoelectron  spectroscopy. 
Especially  important  is  the  variation  of  relative  intensities  as  a func- 
tion of  the  incident  photon  energy  which  makes  possible  the  observation  of 
peaks  absent  for  some  sources,  e.g.,  the  30^  peak  in  acetylene  (Cavell  and 
Allison,  1978).  It  is  obvious  that  no  treatment  of  bound-continuum  tran- 
sitions would  be  complete  without  a concomitant  theoretical  scheme  for  the 
calculation  of  photoionization  cross-sections. 

Extensive  theoretical  studies  of  atomic  photoionization  cross-sections 
are  available  (Bethe  and  Salpeter,  1977;  Marr,  1967;  McGinn,  1970;  Huang 
and  Ellison,  1974a).  Calculations  with  varying  degrees  of  approximation 
for  small  diatomics  (Bates  et  al . , 1953;  Cohen  and  Fano,  1966;  Flannery 
and  Opik,  1965;  Kelly,  1973;  Schneider  and  Berry,  1969;  Tuckwell,  1970) 
and  for  several  hydrocarbons  (Johnson  and  Rice,  1968;  Kaplan  and  Markin, 
1968a,  1968b;  Lohr  and  Robin,  1970)  have  appeared  in  the  literature. 
Calculations  using  uncorrelated  descriptions  of  the  initial  and  final 
states  in  the  CNDO,  MS-Xa  and  HF  approximations  for  the  initial  state  and 
plane  wave  (Kaplan  and  Markin,  1968a,  1968b;  Thiel  and  Schweig,  1971, 

1972;  Rabalais  et  al . , 1974,  1975)  or  orthogonal i zed  plane  wave  (Cavell 
and  Allison,  1978;  Allison  and  Cavell,  1978;  Lohr,  1972;  Rabalais  et  al . , 
1974,  1975;  Huang  and  Ell ison,  1974b)  approximation  to  the  ejected  electron 
have  been  performed.  Erudite  discussion  of  the  theoretical  methodology 
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and  physical  insights  gained  from  these  methods  have  been  presented  by 
several  authors  (Huang  and  Rabalais,  1978;  Lohr,  1972;  Kaplan  and  Markin, 
1968a,  1968b;  Schweig  and  Thiel,  1974;  Ellison,  1974,  1975;  Shirley, 

1978).  Encouraging  results  have  been  obtained  with  correlated  treatments 
of  ionization  cross-sections  using  a plane-wave  approximation  for  the 
ejected  electron  (Purvis  and  Ohrn,  1975;  Herman  et  al . , 1978;  Born  and 
Ohrn,  1979;  Mishra  and  Ohrn,  1980a). 

The  plane  wave  approximation  is  inadequate  for  describing  ionization 
processes  involving  substantial  interaction  between  the  photoelectron  and 
the  ionized  target.  It  also  gives  an  incorrect  angular  distribution  for 
ionization  intensities.  The  uncorrelated  treatment  is,  of  course, 
unsuitable  for  describing  ionization  processes  involving  considerable 
configuration  interaction. 

Recently  many  authors  have  communicated  results  which  show  a complete 
breakdown  of  the  orbital  picture  of  ionization  for  inner  valence  orbitals 
(Schirmer  et  al.,  1977;  Schirmer  and  Cederbaum,  1978b;  Cederbaum  et  al . , 
1978;  Mishra  and  Ohrn,  1980a).  The  sensitivity  of  the  relative  peak 
intensities  to  the  correlation  effects  has  been  pointed  out  by  recent 
publications  (Herman  et  al . , 1978;  Martin  and  Shirley,  1976;  Martin  et 
al.,  1976). 

These  considerations  have  prompted  us  to  put  forth  a correlated 
treatment  using  propagator  amplitudes  to  describe  the  ionization  proces- 
ses and  the  orthogonal ized  plane  wave  (OPW)  approximation  for  the  ejected 
electron.  The  Feynman-Dyson  amplitudes  of  the  electron  propagator 
incorporate  the  correlation  effects  in  ionization  processes  and  lend 
themselves  to  systematic  improvement  (Ohrn,  1976;  Cederbaum  and  Domcke, 
1977;  McCurdy  et  al . , 1977).  The  OPW  approximation  provides  for  the 
necessary  distortion  of  the  plane  wave,  gives  the  correct  form  for  the 


24 


angular  distribution,  and  ensures  the  orthogonality  of  the  initial  and 
the  final  N-particle  states. 

In  the  following  sections  we  present  the  necessary  formalism  followed 
by  the  use  of  calculated  intensities  for  further  investigation  of  the  full 
2p-h  Tamm  Dancoff  Approximation  to  the  electron  propagator. 

2.2  Theory 

The  differential  cross-section  for  bound-free  transitions  induced  by 
incident  radiation  of  angular  frequency  co,  direction  of  incidence  and 
polarization  vector  n,  into  any  number  of  free  states  with  momentum 
p = "hk  in  the  solid  angle  dft  about  t,  using  first  order  time-dependent 
perturbation  theory  and  a semi-classical  description  of  the  radiation-mat- 
ter interaction  is  given  by  (Schiff,  1968;  Smith,  1971). 


da  (£,£ ,n) 
s 01 

dfi 


e2L3k 

27TTTlCa) 


I<N,s|Z 

j 


(2.1) 


where  L is  the  edge  length  of  a very  large  but  finite  cubic  box  used  to 
normalize  the  continuum  functions. 

We  choose  the  many  electron  ground  state  wave  function  of  a closed 
shell  system  as 


l^>  ~ ^0 ( X1 ,x2 ’ ' ' *XN^  (2*2) 

with  x = (r,c)  a compound  space-spin  variable,  and  the  final  state  |N,s> 
is  represented  by 


N,s> 
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Here  and  are  the  two  doublet  spin  components  of  the  (N-l ) - 


Where  is  the  permutation  operator  which  interchanges  electron 
labels  i and  N.  The  form  of  Eq.  2.3  preserves  the  singlet  spin  symmetry 
of  the  system.  Note  that  the  antisymmetric  wave  functions  TQ  and  Ys  are 
exact. 

These  choices  for  |N>  and  |N,s>  lead  to  (Purvis  and  Ohrn,  1975) 


The  first  term  in  Eq.  2.5  relates  the  Feynman-Dyson  amplitude  g (r) 


to  Tq  the  second  term  in  Eq.  2.5  vanishes.  This  is  ensured  by  the  choice 


electron  ion,  v(k,r)  denotes  the  wavefunction  of  the  photoelectron,  and 
the  antisymmetrizer  is 


N-l 


(2.4) 


(2.5) 


where 


N'S  s(r)x(C)  = k(xp.. 


5 • • • 


X|^j_2 » r , l, ) dx^ • • . dx^_2  (2.6) 


with  x(?)  chosen  as  a(c)  or  g(s),  and 


to  the  photoionization  cross-section.  When  v(l<,r)  is  strongly  orthogonal 


v(£,r)  = N0[l  - Il^(r4£)><$£(r4il)|]|£> 


(2.8) 


where 
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|ft>  = L~  ^ exp(ift*r)  (2.9) 

and  ^(r-ft^)  is  the  5-th  CGTO  in  the  computational  basis  set,  centered  at 
position  ft^  in  the  molecular  co-ordinate  system.  For  the  sake  of  clarity 
we  would  like  to  point  out  that  l labels  the  CGTO's,  and  ft.  and  ft  may 
be  the  same  or  different,  depending  on  the  case.  This  orthogonal ization 
unifies  the  OPW  and  the  extended  OPW  formulations. 

It  follows  that  the  normalization  constant  Ng  is  given  by 

Nn2  = {1  - 2L“3  l |<Jt|ft>| 2 + L~3  ££<£|ft><ft|m><&|m>}  1 (2.10) 

l Am 

which  is  taken  to  be  unity  for  a very  large  box. 

These  considerations  reduce  Eq.  2.1  to  the  following  form 
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where 
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Now 


(ft,r)exp(ift  *r)  = L*72[l 
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Retaining  e w , thus,  is  tantamount  to  an  effective  retardation 
of  the  photoelectron.  In  the  following  equations  we  neglect  this  retar- 

~y  -y  -y 

dation,  i.e.,  k-k^  - k (the  dipole  approximation),  whereby 
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(2.14) 


dfl  titticco 


where 


P = v*(l,r)^gs(r)dr 


(2.15) 


The  Feynman-Dyson  amplitudes  gs(r)  are  obtained  from  electron  propa- 
gator calculations  as  a linear  combination  of  some  orbital  basis  (Purvis 


expressed  as  a linear  combination  of  computational  CGTO  basis.  Hence, 
the  amplitudes  may  be  written  as 


The  above  expression  for  the  differential  cross-section,  Eq.  2.14, 
describes  the  ionization  for  a particular  t and  a particular  n in  a 

U) 

specific  photoelectron  direction,  For  comparison  with  spherical  grid 
deflection-type  measurements  (Frost  et  al . , 1968),  done  with  unpolarized 
light  and  a gaseous  sample  of  randomly  oriented  molecules,  we  must 
average  over  all  k^,  n,and  k^.  First  we  present  an  expression  for  total 
cross-section  (a)  suitable  for  comparison  with  spherical  grid  deflection- 
type  measurements  to  be  followed  by  a more  general  formulation  to  obtain 
differential  cross-sections  for  a particular  collection  angle  (e.g.  o ). 

We  first  average  over  all  n and  ^(Doughty,  1964) 


and  Ohrn,  1974)  (e.g.,  the  HF  basis  in  our  work)  which  in  turn  can  be 
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where  n.  and  rL  are  two  polarization  directions  normal  to  k . 

1 L to 
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n2  and  form  a unit  triad  of  vectors. 


IpI2  = 

|nrp|2  + 

|n2*p|2  + Ik^-Pl2 

(2.18) 
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This  integral  may  be  evaluated  by  taking  p as  an  axis  of  a set  of 
polar  co-ordinates  with  k as  the  radial  vector,  and  so 

r O) 
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8tt 


- l^-Pl2)  - fek(l-cos>9)  . \ (2.21) 


which  finally  yields 
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To  evaluate  |p[2  we  define  (Kaijser  and  Smith,  1977) 


5t(t)  = (2irfA 


e~^*r  \(^)dr 


whereby 

P = £[itc*  - mA]exp(-iTc4A)CA(t)L"  /z(2tt)  /z 
and 

Substituting  Eq.  2.24  into  Eq.  2.22  yields 
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where 


Averaging  Eq.  2.26  over  all  photoelectron  direction?  yields 
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which  is  suitable  for  comparison  with  experimental  cross-sections 
measured  in  spherical  grid  deflection-type  spectrometers. 

The  evaluation  of  the  integrals 


exp(-i£4pqHp(£)5*(l<)dfi 


and 


COS0.  . 

mp»k 


exp(-iMpq)5p(lc)cq(lc)dfi 


is  detailed  in  the  Appendix  (A. 21  and  A. 22). 

Our  treatment  of  differential  cross-sections  for  a particular 
collection  angle  is  closely  patterned  after  the  analysis  of  Ellison 
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(1974,  1975).  As  such,  only  skeletal  derivations  are  provided  which  may 
easily  be  fleshed  using  identities  in  the  Appendix.  The  rotation  of  real 
spherical  harmonics  being  cumbersome,  the  treatment  here  is  in  terms  of 
complex  spherical  harmonics  for  notational  compactness.  The  phase 
conventions  are  due  to  Edmonds  (1974)  and  our  treatment  of  the  overlap 
between  the  plane  wave,  and  the  Cartesian  Gaussians  is  in  terms  of 
appropriate  gaussians  in  the  k-space  (Kaijser  and  Smith,  1977)  as  opposed 
to  the  confluent  hypergeometric  functions  (Ellison,  1974,  1975). 

We  define  the  laboratory  co-ordinate  system  (x,y,z)  and  a molecular 
co-ordinate  system  (x'.y'.z1),  their  relative  orientation  being  specified 
by  the  Euler  Angles  = (a,6,y).  Spherical  harmonics  Y£m(e '»<(>' ) of  the 
molecular  co-ordinate  system  are  related  to  those  in  the  lab  system  by 
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and,  conversely, 

VM)  = ( 2 .30) 

a 

where  Dm,m'  are  elements  of  the  usual  rotation  matrix  (Edmonds,  1974). 
Other  relevant  identities  are  tabulated  in  the  Appendix.  The  exciting 
radiation  is  assumed  to  be  directed  along  the  y-axis  with  its  unit  polar- 
ization vector  n along  the  z-axis,  and  the  propagation  vector  for  the 
photoelectron  is  ]<. 

From  Eq.  2.24 

n*P  = Hik  cos0kc^  - Y£]exp(-i(c4£)^£(l<)L"3/2(2Tr)3/2  (2.31) 


where 
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where 
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or,  from  Eq.  2.14, 
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This  expression  provides  differential  cross-section  for  a molecule  with 
its  orientation  fixed  by  the  Euler  angles  B_  = (a,8»Y)  relative  to  the 
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laboratory  co-ordinate  system.  For  randomly  oriented  molecules  we  need 
to  average  Eq.  2.38  over  all  jl.  Furthermore,  for  unpolarized  light,  an 
average  over  all  polarization  angles  is  appropriate. 

The  explicit  dependence  on  B is  displayed  by  the  y's,  ?(t)'s,  and 
exp(-il<*ftpq) , so  we  need  to  treat  the  following: 
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In  the  evaluation  of  Qpq  two  cases  arise.  The  first  case  is  when 
Rpq  = 0.  In  this  case,  using  Appendix  A.  20  and  Eq.  2.29, 
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Q* 

multiplying  by  DQ0(B)  = l and  integrating  over  B,  using  Appendix  A. 15  and 
A. 16, 
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In  the  second  case  Rpq  f 0.  In  this  case,  using  Appendix  A. 2 and  A. 20, 
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Multiplying  Eq.  2.45  by  DQQ  = 1 and  using  Appendix  A. 15  and  A. 16  to 
average  over  all  B_, 
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We  note  that  from  Eq.  2.35  and  Appendix  A. 13  and  A. 14, 
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Also,  from  Eq.  2.35  and  Appendix  A. 20, 
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From  Eq.  2.42  and  2.47, 
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Similarly, 
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In  the  evaluation  of  r two  cases  arise.  The  first  case  is  when 
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Rpq  = 0.  In  this  case,  using  Eq.  2.49  and  Appendix  A. 14,  A. 15  and  A. 16, 
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In  the  second  case,  Rpq  f 0.  Using  Eq.  2.51  and  Appendix  A. 15  and  A. 16, 
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In  the  evaluation  of  x'  two  cases  arise.  The  first  case  is  when 
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Rpq  = 0.  In  this  case,  using  Eq.  2.48  and  Appendix  A. 15  and  A. 16, 


pq 


4tt 


s mimj 


Will.ii 

x YlO(0k’<f)k)  Ylml(0np’cf>np) 
cos0kXpq 


(2.55) 


In  the  second  case,  Rpq  f 0.  Using  Eq.  2.50  and  Appendix  A. 15  and  A. 16, 
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With  these  relationships,  Eq.  2.38  for  differential  ionization 
cross-sections  due  to  polarized  light  with  polarization  vector  n incident 
on  a gaseous  sample  of  randomly  oriented  molecules  may  be  written  as 


da 

dft 


8Tr2e2k 

3mcai 


(I  l [k2cos2ekQ 

p q 


2kcos2ekIm(Xpq)]} 

(2.57) 


where  Tqq  provides  the  spherical  contributions  to  the  distribution  of 
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ejected  electrons.  It  is  well  known  that  the  angular  dependence  of  the 
differential  cross-sections  is  given  by  (Tully  et  al . , 1968;  Yang,  1948; 
Cooper  and  Zare,  1968;  Buckingham  et  al.,  1970;  McGowan  et  al . , 1969) 

H = ??  [1  + BP2(cosek)]  (2.58) 


where  8 is  the  asymmetry  parameter.  Defining 


a 


X 


da 


dfi 


0k=9O° 


and 


da 


dti 


v°° 


from  Eq.  2.59  and  2.60,  we  find 
4-rr(a  +a  ) 

' II  1 

a - 3 


and 


2 (5,,-aJ 

8 = 

(a,+2aA) 


(2.59) 


(2.60) 


(2.61) 


(2.62) 


For  unpolarized  light  the  angle  0^  between  the  polarization  vector 
and  the  propagation  vector  k cannot  be  specified  uniquely.  However,  for 
any  particular  collection  angle,  e.g.  90°  to  the  direction  of  incidence, 
the  angle  between  t and  the  direction  of  incidence  (y-axis)  is  automat- 
ically fixed.  By  defining  a new  polar  co-ordinate  system  with  the 
direction  of  incidence  as  the  polar  axis,  using  Appendix  A. 7,  the  cos0. 

K 


3> 
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dependence  can  easily  be  reformulated.  Finally  averaging  over  all  the 
polarization  angles  gives  (Ellison,  1974,  1975) 

ai  = (2.63) 

which  is  suitable  for  comparison  with  the  cross-sections  obtained  in  the 
most  popular  experimental  setup. 

The  theory  of  molecular  photoionization  cross-sections  lags  far 
behind  the  experiments.  The  present  formulation  meets  the  necessary 
formal  requirements.  It  treats  principal  and  satellite  peak  intensities 
within  a single  theoretical  and  computational  scheme  and  unifies  the  PW, 
OPW  and  the  extended  OPW  approaches  to  the  continuum  wave  functions.  We 
have  presented  two  different  formulations  to  facilitate  an  easy  implemen- 
tation of  this  approach  using  the  two  widely  different  computational 
strategies  (Thiel  and  Schweig,  1971,  1972;  Rabalais  et  al . , 1974,  1975; 
Cavell  and  Allison,  1978;  Allison  and  Cavell,  1978).  The  use  of  single 
determinantal  descriptions  for  ^ and  in  Eq.  2.1  coupled  with  frozen 
orbital  approximation  leads  to  the  OPW  formulation  of  Ellison  (1974, 
1975),  where  the  propagator  amplitude  in  Eq.  2.12  is  replaced  by  the 
corresponding  orbital.  As  such,  this  derivation  also  provides  an  alter- 
nate computational  scheme  within  the  confines  of  HF  approximation. 
Spherical  grid  and  magic  angle  [54.74°  for  which  P2(cos0)  vanishes] 
cross-section  measurements  may  be  used  to  check  the  first  formulation. 

The  measurements  can  ascertain  the  importance  of  the  second 
formulation. 

2 • 3 The  Use  of  Cross-section  Calculations  for  Further  Comparison  of  Full 

2p-h  TDA  and  Other  Self-energy  Approximants 

Real  merit  of  the  scheme  outlined  in  the  previous  section  is  in  its 
unifying  versatility.  The  plane  wave  (PW)  approximation  is  recovered  by 
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retaining  only  first  term  in  Eq.  2.28  and  2.57  (Ellison,  1974,  1975; 

Purvis  and  Ohrn,  1975).  The  correct  angular  distribution 

1(6)  a (1  + |(3cos20-l)}  (2.64) 

given  by  the  OPW  approximation  makes  possible  a study  of  variations  in 
the  angular  asymmetry  parameter  6.  The  parameter  6 depends  on  incident 
ionizing  radiation  and  the  orbital (s)  involved  in  the  ionization  process. 
Accurate  calculation  of  6 thus  serves  as  a potent  assignment  aid  in  the 
interpretation  of  angularly  resolved  photoelectron  spectra.  The  plane 
wave  approximation  leads  to  the  angular  distribution 

1(6)  a cos28  (2.65) 

which  sets  8 = 2 for  ionizations  from  all  orbitals  with  different  angular 
dependence  and  limits  its  utility.  Use  of  propagator  amplitudes  permits 
the  treatment  of  principal  and  satellite  (correlation)  peak  intensities 
within  the  same  scheme.  Summing  up,  the  procedure  outlined  above  could 
make  possible  an  interpretation  and  study  of  finer  details  of  the  photo- 
electron spectrum  in  a way  unactuated  before. 

The  use  of  OPW  approximation  has  definite  advantages.  However,  if 
finer  details  of  the  spectrum  are  not  to  be  elicited  from  the  theory  and 
only  relative  intensities  for  ionizations  by  high  energy  radiation  sources, 
[e.g.,  Mgka  (-1254  ev)]  are  sought,  the  plane  wave  approximation  may  suffice. 

An  inspection  of  Eq.  228  and  2.57  suggests  that  for  high  energy  photoelec- 
trons (large  k),  the  PW  contribution  to  total  intensities  will  be  domi- 
nant, and  the  OPW  corrections  could  be  neglected  since  the  PW  term  is  multi- 
plied by  k3  while  the  OPW  corrections  are  multiplied  by  k or  k2.  The 


Table  V.  Ionization  Cross-sections  and  Potentials  for 
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For  easy  comparison  with  the  Hartree  Fock  and  the  experimental  results,  we  have  normalized  the  relative 
intensities  with  respect  to  the  sum  of  the  cross-sections  for  all  the  satellite  plus  the  principle 
poles  for  the  2 a ionization. 


The  two  deepest  core  and  the  four  highest  virtuals  with  orbital  energies  differing  by  1.0  H and  more  from 
the  valence  region  (-1.70  H - 0.0H)  were  dropped  in  the  solution  of  Eq.  2.29. 

The  two  deepest  core  and  the  three  highest  virtuals  with  orbital  energies  differing  by  1.50  H and  more 
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results  of  Rabalais  et  al . (1974,  1975)  corroborate  this  observation. 

The  PW  approximation  also  has  the  advantage  of  computational  ease. 

The  Feynman-Dyson  amplitudes  determine  the  ionization  intensities. 
Calculated  relative  intensities  using  propagator  amplitudes  may  thus 
serve  as  additional  indicators  of  the  efficacy  of  various  self-energy 
approximants.  In  Tables  V,  VI  and  VII  we  present  the  calculated 
relative  intensities  corresponding  to  the  ionization  of  N2  and  H20  mole- 
cules by  Mgk^  radiation  source,  using  plane  wave  approximation  for  the 
photoelectron. 

An  examination  of  the  relative  intensities,  calculated  in  the 
diagonal  and  the  non-diagonal  approximations,  confirms  the  breakdown  of 
the  orbital  picture  of  ionization  for  the  inner  valence  orbitals; 
however,  no  significant  difference  is  evidenced  for  the  outer  valence 
orbitals.  An  average  of  the  calculated  ionization  potentials  weighted 
by  the  corresponding  cross-sections  gives  the  values  (ev)  32.5,  39.4, 
18.4  for  the  H20  2a^,  N2  2 and  2au,  respectively,  in  close  agreement 
with  the  results  obtained  from  the  diagonal  approximation  and  experi- 
ments. A careful  examination  of  the  results  obtained  from  various 
subsets  reveals  a clear  pattern  in  the  apportioning  of  intensities 
between  the  different  lines,  into  which  the  principal  peak  of  the  orbital 
picture  is  split.  The  use  of  intensities  obtained  from  successively 
larger  subsets  thus  provides  a sound  scheme  for  labelling  the  peaks  of 
the  photoelectron  spectrum. 

The  relative  intensities  obtained  from  various  subsets  differ 
markedly.  Stable  results  are  obtained  only  when  all  orbitals  within 
3.0  H of  the  valence  region  (-1.50, -0.0  H)  are  retained  for  H20  and 
within  2.0  H of  the  valence  region  (-1.70, -0.0  H)  are  included  for  N2< 
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Table  VI.  Ionization  Cross-Sections  and  Potentials  for 
H20  (14  CGTO) 


Orbital 

2fun 

Ediag. 

I.P. (exp. )a 

full 

0 

adiag . 

a(exp. 

2ai 

36.4 

34.3 

32.2 

0.996 

1.000 

1.00 

32.7 

0.004 

3a^ 

14.8 

13.2 

14.7 

0.229 

0.316 

0.260 

lb2 

18.9 

18.2 

18.6 

0.095 

0.109 

0.085 

lbi 

13.1 

11.2 

12.6 

0.137 

0.165 

0.095 

a Siegbahn  et  al . , 1969. 
b tebalais  et  al . , 1974,  1975. 
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The  deepest  core  and  six  highest  virtuals  with  orbital  energies  differing  by  3.0  H and  more  from  the 
valence  region  were  dropped. 
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It  should  be  noted  that  the  sole  use  of  pole  strengths  as  estimates  for 
the  relative  ionization  intensities  in  depicting  a theoretical  photoelec- 
tron spectrum  (Cederbaum  et  al . , 1978)  is  misleading.  For  instance  the 
pole  strength  for  the  liru  pole  of  ^ is  0.92  while  the  relative  intensity 


is  calculated  to  be  0.083. 


CHAPTER  3 


TREATMENT  OF  QUASI-BOUND  CONTINUUM  TRANSITIONS  BY 
DILATION  OF  THE  ELECTRON  PROPAGATOR 


3.1  Resonances 

In  elastic  electron-atom  or  electron-molecule  scattering,  the 
following  process  may  occur: 

e + A — * (A-)*  — *■  e + A (3.1) 

This  involves  a temporary  trapping  of  the  incident  electron  by  the 
target  A,  which  leads  to  severe  distortion  of  the  incident  beam,  causing 
a large  phase  shift  in  the  scattered  wave  which  manifests  itself  as  a 
sharp  increase  of  the  scattering  cross-section.  These  features  character- 
ize resonant  scattering  and  define  resonances  as  quasi -bound  (meta-stable) 
anionic  states  (A  )*  which  decay  by  electron  emission  (Schulz,  1973). 

The  resonance  energy,  i.e.,  the  kinetic  energy  of  the  incident  beam  for 
which  resonances  occur,  is  obtained  from  Eq.  3.1  as 

ER  = K-Ee  = E(A-)*  - EA  (3'2> 

and  is  equal  to  the  electron  affinity  of  the  excited  target  A*.  This 
association  of  resonance  energies  with  electron  affinity  immediately 
suggests  the  use  of  electron  propagator  formalism  for  their  direct 
calculation  in  conjunction  with  a scheme  which  allows  for  the  easy 
differentiation  of  these  quasi-bound  states  from  other  states  in  the 
continuum  part  of  the  spectrum  of  the  A~  Hamiltonian.  Such  a scheme  is 
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presented  by  the  theorems  of  Aguilar,  Balslev  and  Combes  (Aguilar  and 
Combes,  1971;  Balslev  and  Combes,  1971;  Simon,  1972).  These  theorems 
(to  be  referred  as  ABC  Theorems)  provide  a beautiful  characterization  of 
the  pole  structure  of  the  dilatated  resolvent  (Z-Hj-1.  The  dilatation 
transformation  involves  the  replacement  of  r^  by  riexp(i0)  (O<0<^)  in  the 
co-ordinate  representation  of  the  Hamiltonian.  The  Born-Oppenheimer 
molecular  Hamiltonian  is  not  dilatation  analytic  (Simon,  1978),  and  in 
the  following  discussions,  we  limit  ourselves  to  atomic  Hamiltonians. 

The  ordinary  N+l  electron  atomic  Hamiltonian 


1 N+i  , . 

H = 4 l v?  - z y -1-  + Y -L 

2 . “ - l 4-  r . -L . r . . 

i = l i i KJ  ij 


(3.3) 


on  dilatation  transformation  [r^  — *■  r.exp(i0)]  is  transformed  into 


H(0)  = -ie’2i0  l v?  - Ze'i0  l±+  e-ie  l 


i l 


1<j  rij 


(3.4) 


The  ABC  theorems  divide  the  eigenvalue  spectrum  of  H(0)  into  a 
discrete  point  spectrum  ad(0),  a continuous  essential  spectrum  a (0),  and 
a set  of  complex  thresholds  at(0).  The  spectrum  of  H(0)  is  composed  of 
the  following  elements:  a)  A set  of  real  points  in  ad(0)  and  in  at(0). 

The  real  elements  in  ad(0)  are  the  ordinary  bound  states  of  the  (N+l)- 
electron  system;  those  in  a^.(0)  are  the  real  thresholds.  These  real 
elements  in  crd(0)  and  at(0)  are  independent  of  0.  b)  An  essential  spec- 
trum tfg(0)  which  consists  of  a number  of  complex  rays  each  making  an 
angle  -20  to  the  positive  real  axis.  Each  ray  begins  at  one  of  the  real 
thresholds  of  H and  may  be  regarded  as  a rotated  branch  cut.  In 
addition  to  these,  there  also  exists  a set  of  rays  which  begin  at  the 
complex  (resonance)  thresholds  in  at(0).  These  latter  also  make  an  angle 
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-20  to  the  positive  real  axis,  c)  A set  of  complex  discrete  eigenvalues 
in  0^(0)  and  a^.(0).  These  are  defined  to  be  resonances  of  the  (N+l) 
electron  system  with  Z = E - i^-,  where  E and  r are  the  energy  and  width  of 
the  resonance  state,  respectively.  In  the  two-body  case,  the  resonances 
so  defined  have  been  shown  to  be  identical  to  the  poles  of  the  analyt- 
ically continued  scattering  amplitude  (Simon,  1973),  imbuing  this  defini- 
tion with  some  validity.  In  the  general  many-body  case,  the  work  of 
Si  gal  (1978)  is  promising.  Each  resonance  eigenvalue  lies  on  a sector  or 
sheet  of  the  complex  plane  bounded  by  two  of  the  rays  described  in  part  b 
above.  On  a given  sheet,  the  resonance  eigenvalue  is  independent  of  0. 

These  features  make  possible  the  description  of  quasi-bound  continu- 
um transitions  within  the  electron  propagator  formalism.  We  consider  the 
following  eigenvalue  equations  for  the  number-free  dilatated  atomic 
Hamiltonians 


H ( 0 ) 


N> 

0 0 


eV> 

L0'0  0 


H(e)|f>0 


fN±1.N±1 
S,0*S  0 


(3.5) 


where  the  eigenvalue  Eq  for  the  N-particle  ground  state  |0>Q  is  indepen- 
dent of  the  rotation  parameter  0 (Aguilar  and  Combes,  1971;  Balslev  and 
Combes,  1971).  The  complex  conjugate  of  these  exact  eigenfunctions  will 
be  the  eigenfunctions  of  H+(0),  and  the  following  identity  holds  : 


iN±l  = y | N±1  N±1 

£ 1 n 00  n 


(3.6) 


where  the  bar  denotes  complex  conjugation.  In  terms  of  the  dilatation- 
ally  transformed  Hamiltonian  H(0)  e H0,  the  time  development  of 
Heisenberg  operators  is  described  by  the  similarity  transformation 


A(t)  = 


- i Hot 

e 9 A(0)e  9 


A(0)  e A 


(3.7) 
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The  generalized  equation  of  motion  is 

= [A(t),He]_ 


(3.8) 


and  Hg(t)  = H0  is  maintained.  We  define  the  dilated  electron  propagator 
Grs(0,t,t')  = «ar(t)  ; a*(t' )» 

= - ie(t-t')0<N|ar(t)a+(t')i[j>0 

+ i 0 ( t 1 -t ) 0<^| ( t ' )ar(t) |g>Q  (3.9) 


where  0 ( t- t 1 ) is  the  Heaviside  stepfunction. 

In  obvious  analogy  with  the  ordinary  electron  propagator  theory 
(Linderberg  and  Ohrn,  1973),  by  using  identities  Eq.  3.6  and  Eq.  3.7,  the 
Fourier  transform  of  Eq.  3.9  would  provide  the  following  spectral  repre- 
sentation of  the  dilated  electron  propagator 


Grs(e,E) 


«ar’as»E 


"m+  l 

q-K)  n 


,N | | N+l  N+li  +,N 

0 0 1 ar I n ee  n I as 1 0 0 


E-{E®(N+l)-E0(N)}+in 

N. 


<M,  +|N-i>  N-l, 

0 0|as'  n 06  n 1 ar 1 00 


E-{E0(N)-E®(N-l)}-in  J (3*10) 


It  is  obvious  that  the  pole  structure  of  the  dilated  electron  propa- 
gator would  be  similar  to  that  of  the  dilated  resolvents  as  described  by 
the  theorems  of  Aguilar,  Balslev  and  Combes,  except  that  the  poles  here 
correspond  to  energy  differences  between  N and  N±1  electron  states 
(Winkler,  1977;  Donnelly  and  Simons,  1980).  Specifically,  the  dilated 
electron  propagator  will  have  isolated  complex  poles  in  sectors  between 
the  adjacent  branch  cuts,  which  once  uncovered  will  be  independent  of  0. 
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These  complex  poles  z correspond  to 


pN+1 

res  ‘ t0 


rN+l 


z 


res 


(3.11) 


2 


since  the  stable  target  has  zero  width, and  are  to  be  associated  with 
resonances.  It  follows  that,  by  dilating  the  electron  propagator,  a 
computation  performed  on  the  neutral  (N-particle)  parent  would  yield  the 
resonance  energy  and  width  of  interest.  The  construction  of  such  a 
propagator  is  the  main  concern  of  this  chapter.  The  non-hermiticity  of 
dilatated  Hamiltonians  necessitates  some  fundamental  changes  to  the  usual 
scheme  for  the  construction  of  the  electron  propagator.  In  the  next 
section  we  describe  a bi -variational  SCF  procedure  for  dilatated  Hamil- 
tonians. The  use  of  creation  and  annihilation  operators  defined  on  the 
bi -variational ly  obtained  bi -orthonormal  orbitals  in  the  construction  of 
the  dilated  electron  propagator  is  described  in  the  last  section. 

3.2  The  Bi -variational  SCF 

The  typical  dilated  atomic  Hamiltonian  (in  a.u.) 


is  non-Hermitian  for  complex  values  of  the  dilation  parameter  0,  and,  as 
such,  the  variational  theorem  does  not  apply.  However,  there  exists  a 
bi-variational  theorem  (Lowdin,  1978)  for  non-Hermitian  operators.  The 
bi-variational  SCF-equations  for  the  dilated  Hamiltonians  are  derived 
by  extremizing  the  generalized  functional 


ij  i<J  ij 


(3.12) 


(3.13) 
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under  the  constraint  that  the  solutions  $g  and  'Fg  will  be  of  determinan- 
tal  type  and  that  the  constituent  one-electron  orbitals  will  be 
bi -orthogonal . 


$0  = (N!)"^det{(j)1(x1)(j)2(x2)...(|)N(xN)}  (3.14) 

'J'O  = (N ! )-^  det{xp1(x1)^2(x2 ) . . .^N(xN)  > (3.15) 

<+ii*j>  * 5ij  <3-16> 

This  has  been  shown  (Froelich,  1980a)  to  result  in  a canonical  bi-varia- 
tional  SCF  procedure  leading  to  the  effective  Fock  operators 


n 


1(0,4.$)  = ~he'zQ  V*  - e"6^-  + 


1-P 


12 


x2=x2 


r ^ p(x2,x2)dx2 


(3.17) 


and 


where 


+ . 


aT(0,i|»,<J>)  = -%e-20*  V?  - e~6*  -7-  + 

1 rx 


1-P 


x2=x2 


r^  p (x2>x2^dx2 


(3.18) 


P = 


occ 


occ. 


(3.19) 


Extremization  of  the  functional  in  Eq.  3.13  defines  the  necessary 
and  sufficient  condition  (Lowdin,  1981) 


Qp  = pft 


(3.20) 


and  the  following  SCF  equations  result 
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ftrl;.  = e^. 

(3.21) 

(3.22) 

The  N lowest  real  orbital  energies  e ^ determine  the  occupied  spin 
orbitals  which  define  p,  ft,  'Fq,  ft+  and 

The  procedure  followed  uses  a finite  real  basic  set  u = (u^llJLj,  M>N, 
such  that 


i = yc  $ = ud  (3.23) 

(the  bars  on  top  of  ijj  and  $ denote  their  approximate  nature),  and  by 
applying  the  bi-variational  theorem  again  (Lowdin,  1981)  to  the  func- 
tional 


n = <(fr  1 

and  its  complex  conjugate,  we  obtain  the  matrix  equations 


(3.24) 


99  = ACT]  (3.25) 

ft+d  = Adq*  (3.26) 

and  the  bi-orthoganality  relation 

4+Ac  = 1 (3.27) 

where  A = <u|u>  is  the  real  and  symmetric  matrix  and  ft  = <u|ftu>  is 
symmetric  but  complex. 

The  fact  that  H (0)  = H(0*)  = H*(0)  suggests  the  assumption  = 'F* 
and  the  consequent  association  (4>.}  = (Winkler,  1977;  McCurdy  et 
al.,  1980)  which  leads  to  the  functional  of  McCurdy  et  al . (1980).  With 
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these  assumptions  the  following  relations  hold  : 

+ * _+  * 

p = p ft  = ft 

= <u|ft+|u>  = <u* I ft* I u*>  = ft*  (for  real  u)  (3.28) 
1*  t 

This  leads  to  d = c , and  the  modified  bi -orthonormal ity  condition 
c1"  Ac  = 1 holds. 

We  notice  that  if  u were  complex,  then  ft+  f ft*,  d1  f ct,  } = {ip*} 
and  {^. } are  not  bi -orthogonal  any  more. 

The  assumption  {0^.}  = {ip*}  should  be  considered  as  a constraint 
introduced  to  make  the  approximate  many-electron  wavefunction  satisfy  the 
same  relations  as  the  exact  one.  This  choice  permits  the  construction  of 
a bi-variational  SCF  program  driven  by  supermatrices  (Veillard,  1974) 
requiring  a symmetric  density  matrix  and  preserves  the  associated  storage 
and  10  advantages.  It  just  may  be  that  the  extraneous  enforcement  of 
bi -orthogonal ity  condition  is  less  cumbersome  than  the  coupled  iterative 
solution  of  Eq.  3.25  and  3.26  in  a complex  basic  set,  but  it  would  be 
interesting  to  study  the  proximity  of  spaces  { <J>^ } and  where  integral 

driven  SCF  programs  are  available.  A scheme  for  such  a study  is  given  by 

Jankowski  (1976). 

In  our  numerical  investigation  we  have  used  a real  basic  set  of 
symmetry  adapted  CGTO's  and  have  assumed  { cf> ^ } = {^i*}.  With  this  assump- 
tion p = I^qcc^Iqcc  = lyc><uc*|  = |u>cc<u|.  The  convergence  cri- 

terion Sp  = 0 demands  S(cc)  = 0.  As  such,  agreement  of  the  real  part 
of  the  density  elements  from  successive  cycles  within  1.0  x 10-5  was 
enforced  as  a convergence  criterion.  The  stationarity  of  the  total  energy 
was  found  to  be  a less  sensitive  convergence  criterion.  The  generalized 
eigenvalue  Eq.  3.25  was  solved  using  IMSL  routine  EIGZC  (Kaufman,  1975). 


Table  VIII. 


Orbital 


Energies  for  Beryllium  with  (4s/2p)  Basis  Set 


t 


a = 0.0 

a = 0.1 

a = 0.2 

els 

- 4.72163 

(-  4.7251  -0.0096i ) 

(-  4.7341  -0 . 01 68i ) 

e2s 

- 0.30361 

(-  0. 30517-0. 0041i) 

(-  0. 30973-0. 0064i) 

e3s 

0.40697 

( 0.3858  -0.12691) 

( 0.3245  -0.2399i ) 

e4s 

5.48157 

( 5.191  -1.854i  ) 

( 4.333  -3 . 557 i ) 

£2p 

0.05627 

( 0. 05438-0. 0096i) 

( 0. 04879-0. 0178i) 

3p 

0.28766 

( 0.2782  -0.0686i ) 

( 0.2504  -0. 1312i ) 

ETot. 

-14.5683 

(-14.5743  -0 . 00457 i ) 

(-14.5865  -0.0227i ) 

t See  Table  XI  for  the  description  of  this  Basis. 
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Table  VIII  - Extended 


a = 0.3 

(-  4.7437  -0 . 01597 i ) 
(-  0.3171  -0.005i  ) 
( 0.2299  -0.3255i  ) 
( 2.9461  -4.960i  ) 
( 0. 03955-0. 0228i  ) 
( 0.2069  -0.18161  ) 


a = 0.4 

(-  4.746  +0.0064i  ) 

(-  0.3271  +0.002079i ) 
( 0.1130  -0.3725i  ) 

( 1.0945  -0.05911  ) 

( 0.2676  -0.022211  ) 

( 0. 15218-0. 21261  ) 

(-14.5201  -0.16421  ) 


a = 0.5 

(-  2.611  -1 . 526i  ) 

(-  1.627  -8.12191  ) 

( 0. 08388-0. 72921  ) 

( 0.02867-0.12441  ) 

( 0. 05704-0. 100871 ) 

( 0.2694  -0.612751 ) 

(-17.0684  -14.36101 ) 


(-14.5840  -0. 0701i  ) 


Table  IX.  Orbital  Energies  for  Beryllium  with  (10s/6p)  Basis  Set+ 


a = 0.0 

a = 0.1 

a = 0.2 

Is 

. - 4.73323 

(-  4.7335  +0. 000273i 

) 

(-  4.7333  -0. 0000597 i ) 

2s 

- 0.30912 

(-  0. 30913+0. 000042i 

) 

(-  0.30926  +0.0003 06 i ) 

3s 

0.00195 

( 0. 00190-0. 000417i 

) 

( 0. 001764-0. 0008118i) 

4s 

0.00893 

( 0. 00871-0. 00193i 

) 

( 0. 008066-0. 003755i  ) 

9s 

0.8525 

( 0.8242  -0. 2 1 2 5 i 

) 

( 0.7418  -0.4089i  ) 

10s 

6.529 

( 6.210  -2. 119i 

) 

( 5.2682  -0.082i  ) 

2p 

0.00418 

( 0. 00410-0. 0008133i) 

( 0.00387  -0.0016i  ) 

3p 

0.1448 

( 0. 01432-0. 002476i 

) 

( 0.01386  -0. 00496i  ) 

6p 

0.30050 

( 0.2943  -0. 06075i 

) 

( 0.2759  - 0 . 1 1 9 0 i ) 

7p 

0.7898 

( 0.7791  -0. 1737i 

) 

( 0.7207  -0 . 3374 i ) 

'Tot. 

-14.57108 

(-14.5745  -0. 002056i 

) 

(-14.5803  -0.01591i  ) 

t Basis  Set  of  Table  XI  augmented  with  6s-type  GTO's  with  exponents  given 
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Table  IX  - Extended 


a = 0.3  a = 0.4 

(-  4.7293  +0.00022171 ) (-  4.717  +0.00629i  ) ( 

(-  0.3098  +0. 00103i  ) (-  0.3117  +0.00245i  ) ( 

( 0. 001544-0. 0011641  ) ( 0. 001254-0. 001455i ) ( 

( 0. 007036-0. 005376i  ) ( 0. 005685-0. 006708i ) ( 

( 0.6128  -0. 5745i  ) ( 0.4488  -0.6972i  ) ( 

( 3.743  -5.7321  ) ( 1.6988  -0.6920i  ) ( 

( 0. 003493-0. 0023331  ) ( 0. 002976-0. 002986i ) ( 

( 0.01307  -0.007455i  ) ( 0.01179  -0.009963i)  ( 

( 0.2460  -0.17241  ) ( 0.2058  -0.21871  ) ( 

( 0.62742  -0 . 4812 i ) ( 0.5049  -0.5952i  ) ( 

(-14.5719  -0.052821  ) (-14.5125  -0.12031  ) ( 

fo  52l ® 

by  j — I and  4p-type  GTO's  with  exponents  given 
1 n Jn=l 


a = 0.5 

2.5029  -1.54471  ) 

1.0069  -9.42371  ) 

0. 0007396-0. 00154i  ) 
0.003423  -0.006876i) 
0.1724  -0.42211  ) 

0.2213  -1.1055i  ) 

0.002227  -0.003 563 i ) 
0.008219  -0.01349i  ) 
0.2130  -0.39041  ) 

0.5877  -0.125091  ) 

16.2183  -16.79091  ) 

by  (15.56/2.26)"}^. 
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Table  X.  Orbital  Energies  for  Beryllium  with  (14s/16p) 
Basis  Set^ 


a = 0. 0 

a = 

0.1 

a = 

0.2 

£is 

4.73253 

(-  4.73257 

+ 0 . 0000164  i ) 

(-  4.73269  + 

0.0002021  ) 

£2s 

0.30926 

(-  0.30927 

0. 0000107 i) 

(-  0.30928  - 

0 . 00002 1 9 i ) 

*3, 

0.00421 

( 0.00409 

0 . 0009281 i ) 

( 0.00378  - 

0.001801  ) 

0.01957 

{ 0.01904 

0. 004395i  ) 

( 0.01749  - 

0.008531  ) 

^ 13s 

1142.325 

( 1117.799 

- 238.05711  ) 

( 1045.216 

466.2391  ) 

^14s 

7023.166 

( 6877.845 

-1429.885i  ) 

( 6447.731 

2801.6631  ) 

^2p 

0.00007 

( 0.000066 

0 . 00001351  ) 

( 0.000062- 

0.0000265 i ) 

^p 

0.00024 

( 0.000238 

0.0000483i) 

( 0.000224- 

0. 00009471  ) 

^15p 

10.3364 

( 10.1038 

2.17871  ) 

( 9.416 

4 . 261  i ) 

^16p 

25.4107 

( 24.836 

5.4051  ) 

( 23.136 

10.5731  ) 

ETot. 

- 14.57261 

(-  14.57265 

+ 0. 0000697i) 

(-  14.57272  + 

0. 000235i  ) 

t Same  as  that  of  McCurdy  et  al . (1980). 
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Table  X - Extended 


a = 0.3 


a = 0.4 


a = 0.5 


(-  4.7329  < 

(-  0.30931  ■ 

( 0.00327 

( 0.01504 

( 927.528 

( 5750 

( 0.000056 

( 0.000201 

( 8.302 

( 20.385 


0.000826i  ) 
0. 0000616i  ) 
0 . 00257 i ) 
0. 0121i  ) 

675.0681  ) 

40591  ) 

0.0000381  ) 

0.0001371  ) 
6.1541  ) 

• 15.2751  ) 


(-  4.7332  i 

(-  0.3094  • 

( 0.00260  ■ 

( 0.0118  • 

( 769.552 

( 4813.140 
( 0.000047 

( 0.000169 

( 6.809 

( 16.697 


0.002791  ) 

0.0001791  ) 
0.003201  ) 

0.01511  ) 

855 . 843 i ) 
5152.3361  ) 

0.00004891  ) 

0. 000174i  ) 
7.7741  ) 

19i  3 04  i ) 


(-  4.7338  + 0 . 0878i  ) 
(-  0.3098  - 0.0005311  ) 
( 0.00182  - 0.003651  ) 
( 0.00807  - 0.01711  ) 
( 577.62  -10011  ) 
( 3674.476  -6035.401  ) 
( 0.0000368-  0.00005731) 

( 0.000131-  0.0002051  ) 
( 4.999  - 9.04391  ) 
( 12.2338  - 22.47471  ) 


(-14.57260  + 


0. 00064i  ) 


(-  14.5715  + 0.0015731  ) 


(-  14.5672  + 


0.0035631  ) 
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Table  XI.  (4s/2p)  Contracted  Gaussian  Basis  for  Beryllium 


s-type 


Exponents 

Contraction 

Coefficients 

1264.59 

0.001945 

189.37 

0.014835 

43.1591 

0.072091 

12.0987 

0.237154 

3.80632 

0.469199 

1.27289 

0.35652 

3.19646 

0.112649 

0.747814 

0.229506 

0.219966 

1.0000 

0.08231 

1.0000 

p-type 

Exponents 

Contraction 

Coefficients 

6.710 

0.016378 

1.442 

0.091553 

0.4103 

0.341469 

0.1397 

0.685428 

0.0492 


1.0000 
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The  beryllium  atom  was  chosen  as  a first  example  since  the  results  of 
these  calculations  may  be  used  as  input  for  dilated  electron  propagator 
study  of  P shape  resonance  of  Be"  as  electron  attachment  to  2S  Be  ground 
state.  Results  for  basis  sets  I,  II  and  III  (same  as  that  of  McCurdy  et 
al . , 1980)  for  various  values  of  a (in  radians;  a = j-)  are  contained  in 
Tables  VIII,  IX  and  X,  respectively. 

The  eigenvalues  corresponding  to  occupied  orbitals  become  essen- 
tially real  with  enlargement  of  the  basis,  and  the  real  part  of  these 
eigenvalues  is  always  negative,  while  the  eigenvalues  corresponding  to 
the  unoccupied  virtual s are  rotated  into  the  fourth  quadrant  of  the 
complex  plane  irrespective  of  the  basis  set.  The  ground  state  total 
energy,  ETot.  = <$q|H(0) |'Fq>,  becomes  increasingly  real  in  the  same 
fashion.  These  features  suggest  that  the  Hartree-Fock  equations  may  be 
analytically  continued  onto  higher  Riemann  sheets  and  point  to  the 
existence  of  some  self-adjoint  operator  £2'(0)  (=^(0)?)  whose  transforma- 
tion Q(e)  may  be  considered  (Weiner,  1981).  We  note  that  for  smaller 
basis  sets  I and  II  there  is  a complete  breakdown  of  the  dilatation 
analytic  picture  of  fi(9)  for  a = 0.5  radians. 

Surpri s i ngly , the  distinction  between  the  occupied  and  unoccupied 
orbitals  is  maintained  even  with  the  smallest  basis  set  since  the  spec- 
trum remains  strongly  reminiscent  of  the  spectrum  of  rotated  self-adjoint 
operators.  Perfect  bi -orthonormality  was  obtained  in  each  case.  The  use 
of  real  basis  sets  makes  possible  the  use  of  existing  integral  packages 
without  any  modifications,  and  the  flow  of  logic  for  two-electron  inte- 
gral transformation  is  kept  intact.  The  use  of  bi-variational  SCF  in  a 
reasonably  good  basis  set  capped  by  the  construction  of  a dilated 
electron  propagator  in  the  bi -orthonormal  orbital  basis  set  so  generated 
for  the  calculation  of  shape  and  Feshbach  resonances  thus  becomes  a 
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powerful  and  tempting  alternative  to  other  schemes  (McCurdy  et  al . , 1980; 
Donnelly  and  Simons,  1980).  The  construction  of  the  dilated  electron 
propagator  is  described  in  the  next  section. 

3.3  The  Dilated  Electron  Propagator 

The  bi-variational  SCF  furnishes  the  bi -orthonormal  sets  {<(>.}  and 
and  N occupied  orbitals  determine  and  'Fq.  For  the  ideal  case  of 
complete  spin-orbital  bases  {<f>^ } and  (ijj . } the  closure  relation  becomes 

6(x-x‘)  = l (x * )4Jk(x)  (3.29) 

k 

x being  a compound  space  spin  coordinate. 

We  can  write  the  appropriate  zeroth  order  density  operator  as 

p0  = IV<0ol  (3.30) 

with 

Trp0  = 1 (3.31) 


and 

E = <H>  = Tr{HpQ}  (3.32) 

being  stationary.  The  exact  ground  state  solution  can  be  represented  by 
the  density  operator  p = pQ  + (p  - Pq)  (Nehrkorn  et  al . , 1976).  The 
space  of  all  antisymmetric  tensors  of  different  ranks  corresponding  to 
state-vectors  representing  different  numbers  of  particles  generated  by 
the  set  {^},  under  natural  linear  operations  and  exterior  multiplica- 
tions A which  increase  the  rank  of  the  tensors,  constitutes  a Grassmann 
Algebra  G.  We  also  introduce  the  adjoint  algebra  G*,  and  in  terms  of  the 
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duality  <,>  between  G and  G*,  the  interior  multiplication  of  the  left  J, 
and  of  the  right  L,  are  defined  as  follows: 

<x,aiAn>  = <x(aiA)+,n>  = <xja.j,n>  (3.33) 

with  x,  xJa.j  e G*  and  n>  a. An  e G.  Similarly  we  have 

<xAa1,n>  = <x,(Aa1)+n>  = <x,a1Ln>  (3.34) 

x,  xAa1  e G*  and  n,  aYri  e G. 

The  interior  multiplications  a\  and  Ja^  decrease  the  rank  of 
tensors  in  G and  G*,  respectively.  It  is  obvious  that  the  exterior 
multiplications  a^A  and  Aa1  and  the  interior  multiplications  a\  and  Ja^ 
may  be  formalized  as  the  familiar  creation  and  annihilation  operators  on 
the  direct  {i^.}  and  adjoint  spaces,  respectively,  except  that  now 
the  operators  and  their  adjoints  do  not  have  to  have  the  same  domain. 
Specifically  (a^A)+  = Ja^.  These  concepts  are  part  of  the  general 
theory  of  Dual  Grassmann  Algebras  (Froelich,  1980b).  We  cite  some  other 
relevant  results: 

(aiA)+  = jai  (a1 L )+  = Aa1  (3.35) 

[aY,ajL]+  0 = [aiA,a:jA]+  [aiA,aJL]+  = 8..  (3.36) 

We  define  the  electron  field  operators  $(x)  for  the  direct  space 

i(x)  = l ^(x)  a1L  (3.37) 

= l ^ (x)  (Aa1)1, 


(3.38) 
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and  $+(x)  for  the  adjoint  space 


J+(x)  - l *?(x)  (Ja,)  + 
i 

= l <|)*(x)a^ a 


(3.39) 

(3.40) 


and  in  terms  of  the  duality  between  G and  G*,  generalize  the  regular 
second-quantized  representation  of  the  physical  operators,  e.g.,  for  the 
dilated  Hamiltonian  from  Eq.  3.12,  written  as 


H(0)  = l h.(e)+Js£  l g.  -(e) 

i J7i  1J 


i=l 


(3.41) 


H(0)  = 


$+(x1)h1(0)J(x1)dx1  + % it(x1)$t(x2) 


x g12( 0 )^(x2 ^^i )dxidx2 


(3.42) 


l l «fylh(0)lV  arAaSl  + ±l  III 


r s 


r r s s 
s',  s. 


x <4>r4>r.|g(0)ksi|»s.>  arAar,Aas  LasL 


(3.43) 


We  note  that  the  number  operator  for  electrons  will  be 


Noo  * I VaA.  ■ l ", 


op  £ -r 


(3.43) 


and  the  occupation  numbers  are 


Tr{pQarAarL}  e <n  > = 


jl  if  ipr  occupied  in 
(0  otherwise 


(3.44) 


We  partition  the  Hamiltonian  in  Eq.  3.43  as 


H ( 0 ) = Hn(e)  + V(0) 


(3.46) 


where 
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(3.47) 


and 


v(e) 


(3.48) 


with 


<*1^1  II W E <4>i^j|g(0)kk^>  - <4>i<l>J.|g(0)|^k>  (3.49) 

The  linear  space  of  fermionlike  creation  and  annihilation  operators 
introduced  in  the  super-operator  formulation  of  the  propagator  equations 
(Goscinski  and  Lukman,  1970)  is  now  to  be  replaced  by  bi -orthonormal 
operator  spaces 


h = {aiL,aiAaJLakL,...} 
h = {JaiJakJaj.Aa1,...} 


9 • • • 


j>k 


(3.51) 


(3.50) 


with  scalar  product 


(K,|hj)  = Tr{[hj ,ht]+p} 


(3.52) 


We  define  super-operator  Hamiltonian 


H(0)  = Ho(0)  + V(0) 


(3.53) 


and  the  super-operator  identity  I satisfying  the  relations  (Goscinski  and 
Lukman,  1970;  Pickup  and  Goscinski,  1973) 
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H0X  = [X » Hq ( 9 ) ] _ VX  = [X,V(0)]_  (3.54) 

IX  = X V-  X £ h (3.55) 

Analogous  to  the  previous  treatments  (Pickup  and  Goscinski,  1973), 
the  Fourier  transformed  equation  of  motion  of  the  dilated  electron 
propagator  (Eq.  3.10)  may  now  be  written  as 

G(0)  = (a | (El  - H(0) )-1a)  (3.56) 

for  the  electron  propagator,  i.e.,  a matrix  with  elements 


Gij(0)  = (a.  |(EI  - H(9))“1|aJ)  (3.57) 

where  ai  = Ja^ , and  aJ  = aJL. 

Approximation  schemes  can  now  be  developed  exactly  as  in  earlier 
work  (Ohrn,  1976,  1978),  and  instead  of  expression  Eq.  3.46,  we  study 

GO)  = (a|h)(h|(EI  - H(0))h)_1(h|a)  (3.58) 

obtained  by  inner  projection  (Pickup  and  Goscinski,  1973). 

The  zeroth  order  expression  for  the  electron  propagator  is  obtained 
by  using  pg  from  Eq.  3.30  as  the  density  operator, 

h = h1  = (a1*  L}  and  h = h1  e {Ja^}  (3.59) 

in 


G°(0)  = ( a | ( El  - H(0) )a)-1 


(3.60) 


i.e.. 


G°,(e) 


«i3/(E  - ep 


(3.61) 
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This  demonstrates  the  equivalence  of  this  development  to  that  of  the 
unrotated  propagator.  The  structural  equivalence  at  any  order  of  pertur- 
bation theory  between  the  dilated  propagator  6(0)  and  the  unrotated  case, 
demonstrated  by  the  corresponding  choice  of  projection  manifolds  and 
density  operators,  guarantees  the  ease  with  which  the  dilated  propagator 
coalesces  into  the  undilated  one  for  0 = 0 at  any  level  of  approximation. 

For  example,  choosing  p = pQ,  h = h © h , and  h = h © h , with 

h3  = (a  AapLaqL,a  AaaLa3L}  (3.62) 

LXr  U 

(a<6  labeling  occupied,  p<q  labeling  unoccupied  orbitals,  and  i , j label- 
ing unspecified  orbitals)  and 

h3  = {Ja  Ja  Aaa,-laRJa  Aap}  (3.63) 

Hr  P OC 

we  can  write  (Pickup  and  Goscinski,  1973) 

G_1(e)  = El  - e - 1(0)  (3.64) 

where  the  self-energy  matrix  is 


Z(0)  = (a | H(0)h3) (h3 | (El  - H(0))h3)_1(h3|H(0)a) 


(3.65) 


Through  second  order  in  electron  interaction  we  obtain 


Ef  j<e) 


= hlll  ^qMW^AllW 


a p q 


+ hUl 

a 6 P 


E+e  -e  -e„ 
a p q 


E+e  -e  -e0 

p a 3 


(3.66) 


This  expression  for  the  second  order  self-energy  differs  from  that  of 
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Donnelly  and  Simons  (1980)  because  of  a different  partitioning  of  the 
Hamiltonian.  The  expression  of  Winkler  (1977)  was  derived  with  the 
assumption  that  (cp-j ) = and  the  so-called  c-symmetry  (Moiseyev  et 

al . , 1978)  was  employed.  Yet,  the  self-energy  expression  of  Winkler 
does  not  revert  to  the  unrotated  Z?.  for  0=0,  and  could  be  in  error. 

* J 

We  note  that  except  for  the  bi -orthogonal ity,  the  expression  in 

Eq.  3.66  is  remarkably  similar  to  the  unrotated  z? . expression,  and  since 

• «] 

for  0-0,  (4>-j ) = {^-jK  "it  will  be  identical  to  the  standard  expression 
(Pickup  and  Goscinski,  1973).  It  is  obvious  that  retention  of  H instead 
of  its  approximation  by  Hq  in  the  inverse  matrix  in  Eq.  3.65  will  result 
in  the  equivalent  of  the  full  2p-h  TDA  approximation  (Schirmer  and 
Cederbaum,  1978a;  Born  and  Ohrn,  1979;  Mishra  and  Ohrn,  1980a)  for  the 
dilated  propagator.  The  ease  of  extensions  to  mimic  the  higher  order 
decouplings  of  Redmon  et  al . (1975)  or  the  advantages  that  accrue  from 
the  approximation  of  Kurtz  and  Ohrn  (1978)  need  hardly  be  stressed. 

The  approach  outlined  provides  a general  theoretical  format  for  the 
construction  of  the  infinite  order  or  the  renormalized  dilated  electron 
propagator  and  allows  for  simple,  systematic  and  easily  cognizable 
approximation  schemes.  No  assumption  about  the  sufficiency  of  first 
order  correction  to  the  variational  ground  state,  as  a means  of  accommo- 
dating rotation  effects  or  the  applicability  of  the  perturbation  theory 
to  the  perturbation  Hamiltonian  V(0)  = H(0)  - F(0=O),  where  F is  the  Fock 
operator  (Donnelly  and  Simons,  1980)  need  be  made.  The  second  order 
approximation  of  Winkler  (1977)  is  recovered  easily  and  without  invoking 
other  supplementary  concepts  like  the  c-product.  The  scheme  is  compre- 
hensive and  permits  the  treatment  of  cases  where  {4,.}  = is  not  a 

good  approximation  (Froelich,  1980a).  All  the  standard  decoupling 
procedures  for  the  unrotated  propagator  may  be  easily  incorporated. 
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We  notice  that  Eq.  3.64  may  be  rewritten  as 

G"1  (0 ,E)  = El  - L ( 0 , E ) L = e + Z (3.67) 

and  following  Layzer  (1963)  the  poles  of  the  dilated  electron  propagator 
are  the  energy  dependent  eigenvalues  of  the  equation 

L(9,E)Xn(e,E)  = En(e,E)Xn(6,E)  (3.68) 

The  pole  search  scheme  of  Donnelly  and  Simons  (1980)  which  involves  an 
iterative  solution  of  Eq.  3.68  is  quite  satisfactory  and  may  be  incorpo- 
rated in  toto.  The  difference  between  their  scheme  and  the  one  presented 
here  can  now  be  easily  identified  as  the  difference  in  the  choice  of  the 
unperturbed  Hamiltonian  and  the  attendant  orbital  manifold  employed  for 
the  construction  of  L. 

The  recent  results  of  McCurdy  et  al . (1980)  indicate  that  the  shape 
resonances  can  be  treated  at  the  bi-variational  SCF  level.  Feshbach 
resonances  which  involve  more  intricate  correlation  and  relaxation 
effects  are,  however,  poorly  described,  and  the  bi-variational  ASCF  does 
not  furnish  their  lifetime. 

This  success  of  the  bi-variational  SCF  scheme  to  predict  resonances 
points  to  the  optimal  nature  of  the  bi-variationally  obtained  bi-ortho- 
normal orbital  manifold  for  the  construction  of  the  dilated  electron 
propagator.  The  ease  of  incorporating  correlation  effects  to  desired 
levels  would  make  possible  a unified  treatment  of  shape  and  Feshbach 
resonances. 

Our  numerical  investigation  of  the  bi-variational  SCF  procedure 
shows  that  the  CPU  and  storage  requirements  are  comparable  to  those  for 
the  ordinary  SCF  schemes  and  are  entirely  manageable.  Since  the  four 
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index  transformation  packages  (Roos,  1974)  always  utilize  separate 
partial  transformation  of  the  two  indices  on  the  right  followed  by  those 
on  the  left,  the  transformation  of  two-electron  integrals  «X|yv>  over 
the  primitive  computational  basis  set  into  <(p.(p Jip  ip  > would  involve 
minimal  changes  in  the  logic  of  transformations,  and  the  modifications 
to  handle  the  complex  arithmetic  should  not  be  overly  formidable.  Once 
the  transformed  two-electron  integrals  are  available,  due  to  the  ease  of 
the  attendant  order  analysis,  the  construction  of  Z,  and  consequently  L, 
would  be  much  simpler.  It  may  so  happen  that  with  increased  computation- 
al experience  in  the  implementation  of  this  scheme,  one  would  discover  a 
few  energy  dependent  eigenvalues  of  L to  be  entirely  stable  with  respect 
to  changes  in  0>0res,  and  the  need  for  extensive  0-trajectory  calcula- 
tions could  be  supplanted  by  the  ease  of  iterative  solution  of  Eq.  3.68 

for  the  appropriate  symmetry  block  of  L for  a few  0>0  and  identifica- 

res 

tion  of  these  stable  eigenvalues  with  shape  and  Feshbach  resonances.  The 
speed  with  which  iterative  solution  of  Eq.  3.68  was  found  to  converge 
(Donnelly  and  Simons,  1980)  is  most  encouraging. 

Summing  up,  notwithstanding  the  slowness  of  the  complex  arithmetic 
involved,  this  procedure  should  certainly  be  manageable  for  the  kind  of 
systems  to  which  this  methodology  is  being  employed  and  is  warranted  by 
the  fundamental  role  of  bi -orthogonal ity  in  non-hermitean  problems. 

Innate  rigour  and  wider  applicability  should  further  compensate  the  loss 
of  speed. 


EPILOGUE 


The  principal  objective  of  this  dissertation  has  been  the  formulation 
of  theoretical  schemes  for  a correlated  treatment  of  bound-continuum 
transitions.  The  schemes  fabricated  obviate  many  limitations  of  earlier 
attempts  and  could  serve  as  tools  for  an  improved  understanding  of  various 
aspects  of  the  bound-continuum  transitions  within  a general  and  rigorous 
theoretical  framework.  While  the  full  2p-h  TDA  has  been  studied  in  some 
detail  and  criteria  for  its  efficacy  and  dependability  have  been  estab- 
lished, the  OPW  and  the  dilated  electron  propagator  formulations  have  not 
been  investigated  numerically.  The  numerical  groundwork  for  these  two 
schemes  has,  however,  produced  encouraging  results.  A detailed  computa- 
tional investigation  of  these  schemes  is  an  obvious  extension  of  this 
work. 

Throughout,  the  formalism  assumes  closed  shell  reference  states. 

The  open  shell  reference  states  may  be  described  within  the  framework 
established  here,  if  a spin  unrestricted  Hartree  Fock  description  for 
the  reference  state  is  deemed  adequate.  This  less  satisfactory  feature 
may  be  mitigated  by  the  use  of  a unitary  group  approach  to  the  open  shell 
electron  propagator  proposed  recently  (Born  and  Shavitt,  1980).  Due  to 
the  non-dilatation  analyticity  of  the  Born-Oppenheimer  Hamiltonian,  no 
attempt  at  the  description  of  molecular  resonances  has  been  made.  Some 
recent  schemes  (McCurdy  and  Rescigno,  1978;  Simon,  1979)  seem  promising 

in  this  connection,  and  results  obtained  (Rescigno  et  al.,  1980)  are  very 
encouraging. 
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APPENDIX 


We  have  used  the  following  identities  (Harris,  1973) 


00  £ 


exp(±iMpq)  = 4tt  l (±i)A  jo(kRnri)^(7m(eu,d>u)V'Om(0 1 ,4> ' ) 


£=0  m=-£  £ ^ k k ^ 


= 4,1  jo  X {ti)l 


(A.  1) 
(A. 2) 


where  ^£m  anc*  ^£m  c*eno1:e  rea^  anc*  complex  spherical  harmonics  (Edmonds, 
1974),  respectively;  (9k^k)  and  (e 1 ,4>' ) describe  the  directions  oft  and 
^pq’  respectively.  When  t coincides  with  the  molecular  z-axis  (A.l) 
reduces  to 


exp(iiMpq)  = (4^  I (±i)*'j),(kRpq)(2Wl)%yek,lfk) 

£ 

V = Y 
y£0  t£0 

Implicit  in  this  reduction  is  the  addition  theorem 

M“se>  = 2ifrJ_4 

where  9 is  the  angle  between  directions  represented  by  ( 0 1 ,cf> 1 ) and 
(9",<f>").  It  follows  that 


(A. 3) 

(A. 4) 
(A. 5) 
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cose  = P^cose)  = ^ j ^lrn(0, )^1m(0"»4>,,) 


1 


(A. 6) 


fmlj  Yta(6'^')¥im<e"-*")  <A-7) 


The  Clebsch-Gordan  coefficients  are  defined  by  (Harris,  1973) 


m^m^m 


Vi(a)S"i2<n)',««i(a)d 


(A. 8) 


mm^n^ 

'ZX.^2 


yM^^YZ1m^^yi2m2^^d 


(A. 9) 


and  the  orthogonality  relation  for  the  spherical  harmonics  is  given  by 


yo  m (n)^0  m (fi)dfi  = S„  „ 6 

^2m2  2 rn^irir 


(A. 10) 


yl  m Wy!  m ^)dQ.  = 


(A. 11) 


v;m(a)  - (-)mY 


A,-m 


(A. 12) 


The  coupling  rule  for  the  rotation  matrices  is 


„ m1m.:)m1+ml  mlmlm-I-Hii,!,  . 

0 1 : T r w 1 2 r 1 2 1 2 nJ 

m^,m|  n^.m^  j AjA 2 J m1+m2>m{+ni2 


(A. 13) 


and  the  complex  conjugate  of  a matrix  element  is  given  by 


D*! 

mim 


D°*  = i 
u00  1 


(A. 14) 


The  following  integrals  are  useful: 
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sinSdBdcrfy  D**v(B)D*in(B)  - (A.15) 


singdgdady  = 8tt2 


(A. 16) 


The  Fourier  transforms  (Kaijser  and  Smith,  1977)  ££(ft)  may  be 
written  as 


«*<*>  ■ F*<k>  Kwvv 


i n 


= F*<k>  pW  m/VV 


(A. 17) 
(A. 18) 


whereby 


F,WF>)  (ek,^.mAA) 

J 11  J J 

o m 1 J Mm.m. 

Fi(k>  1 1 cicj  ^1^  ! WW 


(A. 19) 


= F,Ck)Fm{k)  | I c,c-  V- (ek.*k)Y.  (ek,*k) 

' J ' ' J J 

. i J n Mm.m. 

FZ^Fm^  n,t.  . , cicj  CL£.£. 


i j L-I^-AJ  1 J LVj 


(A. 20) 


We  evaluate  the  integral  exp(-i£*ft  )e  (k)£*(k)dft  for  a linear 

rH  r H 

molecule  with  ftpq  along  the  z-axis.  For  non-linear  molecules,  similar 
expressions  apply  employing  appropriate  rotation  matrices  (Schweig  and 
Thiel,  1974).  The  integral  is  easily  evaluated  when  ft  = ft,  so  we  treat 
the  case  ftpq  f ft  using  (A. 3),  (A. 8),  (A. 17)  and  (A. 19)  : 
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exp(-iMpq)?p(lc)Cq(t)dfl 


£ .+£ . 

(4u)‘5  F (k)F*(k)  l l l i 

P q * J HV'j 


n n p i om.m. 

* (A.2i) 

■ J 


We  evaluate  the  integral  cos0+  exp( -i£4nJ  £(£)£*(£)  dft  using 

J rn  |<  PM  P q 

(A. 3),  (A. 6),  A. 8),  (A. 17)  and  (A.19): 


COS0fn  ,t  exP(_i^’^pq)^p(  k)5q(k)dn 
P 


£.+£. 
1 J 


- ^ Fp(k)  n . . i 


i J L=|£i-£. 


L+l 


l , (-1)£  j0(kR)(2£+l) 


£=|L-1| 


£'  pq' 


where 


nm.m . 


y y (0  d)  ) rnnO  r 1 J 

n=-l  ln  y V LU  Ui‘1j 


(A. 22) 


0-».  = cos 


-1 


mp 

l*pxl 


(A. 23) 


m*  ' 


iju.  = tan-1 
rfL 


-*x 

mpj 


(A. 24) 


The  superscripts  x,y,z  denote  the  various 


components  of  the  vector 
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